January 28, 2013 



ICCUB-11-145 
MAD-TH- 11-03 



Thermodynamics and Instabilities 

of a 

Strongly Coupled Anisotropic Plasma 

David Mateos 1,2 and Diego Trancanelli 3,4 

1 Institucio Catalana de Recerca i Estudis Avangats (ICREA), Passeig Lluis Companys 23, 

E-08010, Barcelona, Spain 

2 Departament de Fisica Fonamental (FFN) & Institut de Ciencies del Cosmos (ICC), 
Universitat de Barcelona (UB), Marti i Franques 1, E-08028 Barcelona, Spain 

3 Department of Physics, University of Wisconsin, Madison, WI 53706, USA 

4 Department of Physics, University of California, Santa Barbara, CA 93106, USA 

dmateosOicrea. cat , dtrancanOphysics . wise . edu 

Abstract 

We extend our analysis of a IIB supergravity solution dual to a spatially anisotropic finite- 
temperature N = 4 super Yang-Mills plasma. The solution is static, possesses an anisotropic 
horizon, and is completely regular. The full geometry can be viewed as a renormalization 
group flow from an AdS geometry in the ultraviolet to a Lifshitz-like geometry in the infrared. 
The anisotropy can be equivalently understood as resulting from a position-dependent 6*-term 
or from a non-zero number density of dissolved D7-branes. The holographic stress tensor 
is conserved and anisotropic. The presence of a conformal anomaly plays an important 
role in the thermodynamics. The phase diagram exhibits homogeneous and inhomogeneous 
(i.e. mixed) phases. In some regions the homogeneous phase displays instabilities reminiscent 
of those of weakly coupled plasmas. We comment on similarities with QCD at finite baryon 
density and with the phenomenon of cavitation. 
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1 Introduction 

A remarkable conclusion from the Relativistic Heavy Ion Collider (RHIC) experiments [1, 2] 
is that the quark-gluon plasma (QGP) does not behave as a weakly coupled gas of quarks 
and gluons, but rather as a strongly coupled fluid [3, 4]. This renders perturbative methods 
inapplicable in general. The lattice formulation of Quantum Chromodynamics (QCD) is also 
of limited utility, since for example it is not well suited for studying real-time phenomena. 
This has provided a strong motivation for understanding the dynamics of strongly coupled 
non-Abelian plasmas through the gauge/string duality [5, 6, 7] (see [8] for a recent review of 
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applications to the QGP). The simplest example of the duality is the equivalence between 
four-dimensional N = 4 super Yang-Mills (SYM) theory with gauge group SU(N C ) and type 
IIB string theory on AdS§ x S 5 . We recently extended this duality to the case in which the 
SYM plasma is spatially anisotropic [9]. 1 In this paper we provide additional details and 
extend our analysis of the anisotropic supergravity solution. Previous holographic studies 
of anisotropic plasmas include solutions dual to non-commutative SYM theories at finite 
temperature [10, 11, 12, 13, 14], gravity duals of p-wave superconductors[15, 16, 17, 18, 
19, 20, 21], solutions dual to M = 4 SYM theory coupled to electromagnetic fields with 
[22, 23, 24, 25] or without [26, 27] flavour, and anisotropic solutions with singular horizons 
[28]. 

Part of the motivation for our work comes from the fact that the plasma created in a 
heavy ion collision is anisotropic. The system is locally anisotropic for some short time after 
the collision, r < r iso , and becomes locally isotropic afterwards. An intrinsically anisotropic 
hydrodynamical description has been proposed to describe the early stage after the collision 
[29, 30, 31, 32, 33, 34, 35, 36, 37], in which the plasma is assumed to have unequal pressures 
in the longitudinal and transverse directions. We will encounter this feature too. At a time 
r iso the pressures become approximately equal and a standard hydrodynamic description can 
be applied in which the stress tensor is assumed to be locally isotropic, meaning that each 
little cube of QGP is isotropic in its own rest frame. Even in this phase, certain observables 
may be sensitive to the physics in several adjacent cubes, in which case there is no frame in 
which they can be computed as if no anisotropy were present. A simple example is the drag 
experienced by a heavy quark moving through the plasma. On the gravity side the quark 
is modeled by an infinitely long string [38, 39], which is therefore sensitive to the physics in 
a large region of the plasma. On the gauge theory side the size of the string corresponds 
roughly to the size of the gluon cloud that dresses the quark. 

Further motivation for our work is provided by the fact that weakly coupled plasmas, both 
Abelian and non-Abelian, are known to generically suffer from instabilities in the presence 
of anisotropics [40, 41, 42, 43, 44]. In fact, it has been proposed [45, 46, 47, 48, 49, 50, 51] 
that these instabilities might be responsible for the seemingly small value of r iso < 1 fm. It is 
therefore interesting to understand whether this type of instabilities also occur in anisotropic 
plasmas at strong coupling. We will see that the gravity solution that we will present exhibits 
instabilities reminiscent of weak-coupling instabilities. 

Another interesting aspect of our solution also related to the QGP is the fact that it is 
sourced by a position-dependent #-term — see eqn. (1) below. It has been suggested that 
'bubbles' or regions with a spacetime-dependent theta-angle, 9(t, x), may be formed in heavy 
ion collisions [52, 53, 54, 55, 56, 57]. In QCD an effective 0-angle of this form is equivalent 
to a non-zero axial vector potential = 8^9/ 2N t that couples to the axial current as 
■tp^f^^TpA^. This leads to interesting effects, including the generation of an electric current 
along a magnetic field — the so-called chiral magnetic effect [52, 53, 56, 57]. 

We emphasize that the purpose of this paper is to study in detail the anisotropic SYM 

1 The term 'anisotropic' is sometimes used in the literature to refer to systems in which there is an asymmetry 
between time and space, i.e. systems with no Lorentz invariance. In this paper we will use the term 'anisotropic' 
to refer exclusively to spatially anisotropic systems, i.e. to systems in which there is an asymmetry between 
different spatial directions. 
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plasma at strong coupling through its dual gravity description. Although we will comment 
on superficial similarities and differences with the physics of the QGP, an in-depth analysis 
of the extent (if any) to which a close connection with the QGP exists is beyond the scope of 
this paper. For example, we leave for future work questions such as whether the instabilities 
we will identify can be thought of as a strong-coupling counterpart of the weak-coupling 
instabilities mentioned above. 

At a more theoretical level, motivation for our work is provided by the fluid/gravity cor- 
respondence [58, 59, 60, 61, 62, 63, 64, 65] (see [66, 67, 68] for reviews), which asserts that 
the long-wavelength effective dynamics of black brane horizons in asymptotically anti-de 
Sitter spacetimes is equivalent to the long-wavelength effective dynamics of the dual finite- 
temperature plasma, namely to fluid dynamics. It would certainly be interesting to extend 
this correspondence to the case in which the fluid is allowed to be intrinsically anisotropic, 
for which purpose our solution may provide a starting point. 

An interesting connection also exists between our work and the so-called blackfold approach 
[69, 70, 71, 72] to black hole dynamics, of which the fluid/gravity correspondence may be 
viewed as a holographic implementation. In this framework the long- wavelength dynamics of a 
black hole horizon is described by a hydrodynamic theory living on a dynamical worldvolume. 
The black hole is embedded in an arbitrary spacetime, which in particular is not necessarily 
asymptotically AdS, and may carry dipole charges which can be thought of as associated with 
p-dimensional objects dissolved in the black hole's horizon. In these cases the corresponding 
effective dynamics describes an anisotropic fluid with a conserved p-brane charge [72] . As we 
will see, the anisotropy in the SYM plasma that we will study can be thought of as resulting 
from a dissolved 2-brane charge. 

A final motivation for our work comes from the possible applications to condensed matter 
systems. In particular, our solution may be viewed as an ultraviolet (UV) completion of a 
charged, anisotropic Lifshitz-like infrared (IR) geometry, as already observed by the authors 
of Ref. [73] . This allows for a rigorous holographic study of the Lifshitz geometry and provides 
a microscopic description of the dual theory. 2 More generally, anisotropic strongly coupled 
condensed matter systems exist in Nature (e.g. liquid crystals) and we hope that our work 
may contribute to the application of holographic methods to this type of systems. 

The type IIB supergravity solution that we will present is a finite-temperature generali- 
zation of that of Ref. [73], which provided much inspiration for our work. In the Einstein 
frame the ten-dimensional geometry factorizes as A4 x X, where X is a five-dimensional Ein- 
stein manifold which we will take to be X = S 5 for simplicity, and Ai satisfies the following 
properties: 

1. It is static and anisotropic (along the gauge theory directions). 

2. It possesses a horizon and it is regular on and outside the horizon. 

3. It approaches AdSs asymptotically. 

The staticity requirement is motivated by the desire to work in the simplest possible setup. 

2 Embeddings of Lifshitz geometries in string theory were constructed in [74, 75, 76, 77, 78, 79, 80, 81, 82]. 
Recent studies of AdS completions of charged, spatially isotropic Lifshitz geometries include [83, 84, 85, 86, 
87, 88, 89, 90]. 
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In other words, we wish to be able to study the thermodynamics of the system, its response 
when it is slightly perturbed, etc. The presence of a horizon is dual to the existence of a 
finite-temperature plasma in the gauge theory. The requirement of regularity guarantees that 
calculations are unambiguous and well defined. Finally, the AdS boundary conditions ensure 
that holography is on its firmest footing. More specifically, the fact that our configuration 
solves the type IIB supergravity equations of motion and asymptotically approaches AdS§ x S 5 
implies that it is dual to N = 4 SYM theory deformed by a (marginally) relevant operator, 
and thus that it is solidly embedded in string theory. 

As in [73], the deformation of the N = 4 theory that we consider corresponds to the 
addition of a ^-parameter that depends linearly on one of the three spatial coordinates, 
6 = 2irn D7 z, where {t, x, y, z] are the gauge theory coordinates and constant with 

dimensions of energy which (as we will review) can be interpreted as a density of D7-branes 
distributed along the z-direction. In other words, the total gauge theory action takes the 
form 

1 

87T 2 



= SV=4 + SS, SS=^ I 6{z) TvFAF. (1) 



This clearly breaks isotropy (and CP) but not translation invariance, since integration by 
parts yields 

SS oc -n m j dz A Tr (^A A F + ^A 3 ^j . (2) 

Incidentally, this expression shows that, if the z-direction is compactified on a circle, the 
resulting three-dimensional theory contains a Chern-Simons term. If in addition antiperiodic 
boundary conditions are imposed for the fermions around the circle, then the theory flows at 
long distances to a Chern-Simons theory [91]. The deformation (1) breaks all the supersym- 
metries of the four-dimensional SYM theory. Supersymmetry-preserving deformations with a 
space-dependent 0-angle (and coupling constant) have been considered in [92], motivated by 
the construction of Chern-Simons theories with N = 4 supersymmetry in three dimensions. 

The complexified coupling constant of the M = 4 theory is related to the axion-dilaton of 
type IIB supergravity through 

r=— + — =X + ^-K 3 

Thus we expect that the gravity solution dual to the deformation (1) will have a position- 
dependent axion field of the form \ = az where, as we will see, the constant a is given 
by 

An D7 

a= 4^VV (4) 

where A = g^ M N c is the gauge theory 't Hooft coupling. 3 Since the axion is magnetically 
sourced by D7-branes, this suggests that it should be possible to interpret the solution in 
terms of a number of D7-branes dissolved in the geometry [73] — see Fig. 1. We will elaborate 
on this in Section 4. Here we just anticipate that the D7-branes are wrapped on the S 5 factor 
of the geometry, extend along the xy-directions and are homogeneously distributed along the 



3 IIB supergravity solutions with a running axion-dilaton were constructed in [93, 94, 95]. Lifshitz-like 
solutions with a linear axion have been considered in [79, 81, 82]. 
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Figure 1: D7-branes dissolved in the geometry. 



z-direction with uniform density n D7 = dN DT /dz. In view of the relation (4) we will use n D7 
and a interchangeably to refer to the D7-brane number or charge density. The orientation 
of the D3- and D7-branes that give rise to our solution can be summarized in a standard 
notation by the array 
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S 5 




D3 
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X 


X 


X 






n D7 


D7 


X 


X 


X 






X 



where u is the holographic radial coordinate in AdSs. We stress that, since our solution 
incorporates their full back-reaction, the D7-branes are completely 'dissolved' in the geometry, 
just like the iV c D3-branes that give rise to AdS§ x S 5 . We also emphasize that, unlike the case 
of D7-branes used to introduce flavour (quark) degrees of freedom in ftf = 4 SYM [96, 97, 98], 
the D7-branes considered here do not extend in the radial direction. Consequently, they do 
not reach the AdS boundary and they do not add new degrees of freedom to the SYM theory. 

We will verify that translation invariance is indeed preserved by showing that the holo- 
graphic stress tensor is conserved. The calculation of the stress tensor will require that we 
renormalize the gravity action by adding appropriate counterterms, which have been deter- 
mined in [99, 100]. In the process of renormalization we will discover that the deformation 
(1) induces a non-zero conformal anomaly, (T?) ^ 0. This will play an important role in the 
thermodynamics of the system, most crucially because it will imply the existence of an arbi- 
trary reference scale /x, a remnant of the renormalization process much like the subtraction 
point in QCD. This in turn means that the physics does not only depend on a/T, but on 
two independent dimensionless ratios which may be taken to be a/fi and T/fi. A detailed 
study of the thermodynamics of asymptotically locally AdS spacetimes, in particular in the 
presence of a holographic conformal anomaly [101, 102], was presented in [103]. 

As mentioned above, in the UV our solution approaches AdS$ x S 5 , as appropriate for a 
marginally relevant deformation of the N = 4 theory. In this limit the solution is of course 
invariant under the rescaling x 1 — > kx l ,u — > ku. We will see that in the IR it approaches a 
Lifshitz-like solution whose metric is invariant under a scaling as above for all the coordinates 
except for z, which scales as z — > k 2 ^ 3 z. In this sense the solution may be viewed as a 
renormalization group (RG) flow between an isotropic UV fixed point and an anisotropic IR 
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fixed point. 4 The zero-temperature version of this flow was found in [73], and in this case the 
string-frame metric exhibits a naked curvature singularity deep in the IR (see Section 6 and 
Appendix B for details). However, we emphasize that at finite temperature the singularity is 
hidden behind the horizon and the solution is completely regular on and outside the horizon, 
exhibiting no pathologies regardless of the number of D7-branes, which in particular can be 
larger than 24. Incidentally, this provides a completely explicit counterexample (at least in 
the presence of a negative cosmological constant) to the somewhat extended belief that 24 
is an upper bound on the possible number of D7-branes. We will comment further on this 
point below. The finite-temperature deformation of the IR Lifshitz geometry was also found 
analytically by the authors of [73] , which allowed them to show that in this case the entropy 
density scales as s ~ a 1 / 3 ?" 8 / 3 . Holographic renormalization for supergravity solutions dual 
to field theory RG flows was initiated in [104, 105]. 

The organization of the paper is as follows. In Section 2 we review the solution presented in 
[9] and discuss the entropy density. In Section 3 we compute the stress tensor particularizing 
the results of [99, 100], and discuss the role of the conformal anomaly. In Section 4 we discuss 
the definition of the chemical potential conjugate to the D7-brane number density n D7 and 
review several aspects of the thermodynamics of anisotropic systems. In Section 5 we apply 
the results from previous sections to the construction of the phase diagram in the (o, T)-plane 
(i.e. in the canonical ensemble). We identify regions in which the homogeneous anisotropic 
phase is stable, metastable or unstable. In the latter two cases the stable configuration is 
a mixed (inhomogeneous) configuration in which the anisotropic and the isotropic phases 
coexist. A summary of our results is contained in Section 6, where we also elaborate on 
several topics including the similarities and differences of our system with weakly coupled 
plasmas and with QCD at finite baryon density. 

We conclude the paper with a series of appendices in which we collect some of the more 
technical details of our computations. In Appendix A we write down the equations of motion 
and describe the procedure we have used to solve them. In Appendix B we compare our 
parametrization of the geometry with the parametrization employed in [73]. In Appendix C 
we derive the relations between thermodynamic potentials and pressures in an anisotropic 
system. In Appendices D and E we obtain analytic expressions for the metric and the various 
thermodynamical quantities in the regimes T>a and T a, respectively. 

2 Action and solution 

In this section we collect the main results about the solution of interest. Additional details 
can be found in Appendix A. 

The ten-dimensional solution we seek is a direct product of the form Al x S 5 , where the 
radius L of the S 5 is constant in the Einstein frame and is given by 

L 4 = 4irg s NJ 4 s = X£ 4 S , (6) 

where A = <7ym-^c is the gauge theory 't Hooft coupling. Therefore M. can be viewed as a 
solution of five-dimensional supergravity with a negative cosmological constant A = —6/L 2 . 

4 This is a slight abuse of language because, although the metric in the IR is scale-invariant, the dilaton 
and the Ramond-Ramond forms are not. See Appendix B. 
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Since only the metric g^ v , the axion \, and the dilaton will be excited, it suffices to consider 
the five-dimensional axion-dilaton-gravity action, which in the Einstein frame takes the form 



1 

2^ 



M 



g[R + \2 



1 



(d<f>f 



1 



,2<i 



(dx)'' 



+ 



1 

2^ 



- 7 2K, 



dM 



(7) 



where 2k 2 = 16irG is the five- dimensional gravitational coupling. For convenience we will 
measure all lengths in units of L, i.e. we have set L = 1 in eqn. (7), which implies that 
G = 7r/2iV c 2 and k 2 = Air 2 /N 2 (see e.g. [8]). The solution takes the form 



ds 2 



X 



e 2< 



az . 



-TB dt 2 + dx 2 + dy 2 + Udz 2 + 

4> = 4>(u) , 



du 



(8) 
(9) 



where J 7 , B and T~L are functions of the holographic radial coordinate u. Reparametrization 
invariance has been used to fix g xx and g yy in (8); once this is done, B cannot be set to 
unity in general. The solution is clearly isotropic in the xy-directions, but not in the z- 
direction, unless % = 1. The axion, which is dual to the gauge theory #-term, is responsible 
for inducing the anisotropy. J 7 is a 'blackening factor' that will vanish at the position of the 
horizon, u = ii H . The boundary of the spacetime is at u = 0. 

The functions J 7 , B and % only depend on the radial coordinate u and are given in terms 
of the dilaton as 



n 

T 

& 
B 



e 2< 



b' + utf') 
1 



a 2 ei^(4u + u 2 <j)') + 160' 



7 (240' - 9u0' 2 + 2Ou0") 



(10) 
(11) 
(12) 



24 + Wu<p' 

where the primes denote differentiation with respect to u. Eqn. (12) determines B up to an 
arbitrary multiplicative constant. Eqns. (8)-(12) are a solution of the field equations provided 
the dilaton obeys the third-order equation 







2560V" - 160' 3 (7ucf)' + 32) 



7,:, 



+ 



ua 2 e^ (ucjJ + 4) + 160' 



u (hup + 12) (u(f>" + <p') 



x 



13uy 4 + 8-u(lluV' 2 



+2u<p 



600" - 
5«y" + 53n 2 0" + 36) + 



12110"') + u 
6' f30uV' 2 



2 j/3 



64mV" 



>" + 96) 



288 + 32u 2 



.(13) 



Without loss of generality (see Section 4.1), we will seek solutions such that the dilaton 
vanishes at the boundary, i.e. such that 4> bdTy = 0(0) = 0. Eqs. (10) and (11) expanded near 
the boundary then imply ,F bdry = H hd r y = 1- We will also choose the arbitrary constant in 
B such that B hdiy = 1. With these normalizations the temperature and the entropy density 
are readily computed. We first note that the Euclidean continuation of the metric (8) in the 
(£ E ,u)-plane near u n reads 



dsi 



TiB H (u H - u) dt 2 + 



du 2 



Fi{u n - u) 



(14) 
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where T\ = The standard requirement that the metric be regular at u = u H 

then determines the period of the Euclidean time, 5t E , which we identify with the inverse 
temperature: 

r i (15) 

St E Air v ; 

The entropy density is simply obtained from the area of the horizon. The area element on a 
t = const., u = u H hyper surf ace is 

g-f^H 

dA H = 5 — dx dy dz , 

so the entropy density per unit volume in the xyz-directions is 

-i^ry^kr <16) 

The isotropic black D3-brane solution is a solution of the equations above with a = and 

u 4 1 

= 0, B = n = l, F=l- — , u H = —. (17) 
Eqn. (16) then yields the familiar expression for the entropy density of J\f = 4 SYM [106]: 

s\T) = ^N 2 c T\ (18) 

At this point we can perform an interesting check on our solution. As mentioned above, the 
zero-temperature solution is a domain-wall-like solution in the radial direction, interpolating 
between an AdS geometry in the UV and a Lifshitz-like geometry in the IR [73]. The radial 
position at which the transition takes place is set by the anisotropic scale, a. Thus we expect 
that in the limit T S> a the entropy density should scale as in (18), since in this limit the 
horizon should lie in the asymptotic region where the geometry is approximately AdS. In the 
opposite limit we expect the entropy density to scale as 



tiYc V/3 T 8/3 (lg) 



with c cnt a numerical coefficient, since this is the appropriate scaling [73] in the IR Lifshitz-like 
region in which the horizon lies when T -C a (see Appendix B for details). These scalings 
are precisely reproduced by the entropy density computed with our numerical solution (see 
Appendix A for details), as shown in Fig. 2. This plot was produced by evaluating the entropy 
density for many different values of a and T. We see that for a <C T the points are aligned 
along the horizontal axis, thus reproducing (18). In the opposite regime a>T, the points 
are aligned instead along a straight line with slope 1/3, which means that the entropy density 
scales in this case as in (19). In between, the entropy density smoothly interpolates between 
the two limiting behaviours. 
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Figure 2: Log-log plot of the entropy density as a function of a/T, with s° defined as in eqn. (18). 
The dashed blue line is a straight line with slope 1/3. 



3 Stress tensor 

The energy density and the pressures of the deformed J\f = 4 theory can be obtained by 
calculating the holographic stress tensor, defined as the variation of the supergravity action 
with respect to the boundary metric. The action suffers from large-volume divergences, 
which can be regularized and subtracted by a procedure called holographic renormalization 
(see e.g. [107] and references therein). The total on-shell action takes the form 

Son-shell = Sbulk — Set 5 (20) 

where the counterterm action can be obtained by particularizing the results of Refs. [99, 100] 
to the case of interest here. Since we will be interested in the thermodynamics of the system it 
is convenient to work here in Euclidean signature. In this case S bulk is given by the Euclidean 
continuation of (7) and, omitting terms that vanish identically for the case at hand, the 
counterterm action takes the form 

S ^ = hj d4x ^ ( 3 " \ e ** d #?x) ~ l °Z v J d * x ^ A + |( c -n - 1) J d 4 x^A , (21) 

where A^ij, (f>, x) is the conformal anomaly in the axion-dilaton-gravity system [99, 100]. 
The coordinate v is the standard Fefferman-Graham (FG) coordinate, in terms of which the 
metric near the boundary takes the form 

dv 2 

ds 2 = — s- + jij(x, v) dx l dx j . (22) 



Note that if we were to restore the AdS radius L then the coefficient of the second term in 
eqn. (21) would become \og{v/L). 
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The counterterm action (21) is understood to be evaluated on a constant-f slice, the bulk 
action is understood to be evaluated by integrating from the horizon down to the lower cut- 
off v, and the limit v — > is understood to be taken at the end of the calculation. The 
logarithmic term in (21) is necessary to cancel certain divergences from the bulk action, 
and it breaks diffeomorphism invariance in the bulk. This is the origin of the conformal 
anomaly in the boundary. The third term in (21) is finite and can therefore be added with an 
arbitrary coefficient, which we have written as (c sch — l)/4 to simplify subsequent equations. 
The freedom to add this term is part of the general freedom in the choice of renormalization 
scheme. 

The fields have the following general expansions near the boundary: 

l r 



7y = ~2 



$ + + v A (gff + 2 log v $>) + (v°)\ , (23) 

+ v 2 (j) {2) + v 4 (V 4) + 2 log v (4) ) + 0{v e ) , (24) 
X = X {0) +v 2 x^+v 4 ( X (4) + 2logvx {4) )+O(v 6 ). (25) 



In our case, = r\ij is the flat gauge theory metric. The axion is simply 

x = x (°) = az, (26) 

and the dilaton takes the form 

, a 2 2 (2Ba 47 a 4 \ 4 a 4 



v 2 +[ — - )v 4 v 4 logv + 0(v b ) . (27) 

V 7 4032 / 6 5 v ' y 1 

The expansion of the metric can be obtained from the expansion of the functions J-, B and 

n, 

Ua 2 „ / lla 4 \ 4 7a 4 4l n , fiN 

T = 1 + ~^v 2 + ( ^4 + W J v 4 + ~^v 4 log v + 0(v 6 ) , 

B=l — ^v 2 +(b 4 - v 4 - —v 4 \ogv + 0(v G ) , 



24 V 144 / 12 

n = 1 + a \v 2 - ( 2 -f - ^ j v 4 + ~-v 4 log, + 0(„") , (28) 



and it reads 



1 a 2 2 ( 23S 4 2,F 4 2749a 4 a 4 , 4 

to = -H f + H 1 2log?; W; 4 + 0('(; b ), 

y 24 \ 28 4 16128 48 

a 2 9 / 5^4 Fa 71a 4 a 4 , \ , , ; 
Q w = 1 v 2 + i - + 2bgvU 4 + Ot) 6 ^ 



24 V 28 4 1792 48 
5a 2 2 / 13^4 Ta 1163a 4 



y 24 V 28 4 16128 16 J V 7 V 7 

^4 (a, T) and 64(0, T) are two coefficients that are not determined by the asymptotic equations 
of motion, as expected, but that can be read off once a particular solution is known. Below we 
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Figure 3: Expansion coefBcients Ta and B4. In (a) we fix T = 1 and plot J-4/T 4 and K4/T 4 as 
functions of a/T. In (b) we fix a = 1 and plot J^/a 4 and £? 4 /a 4 as functions of T/a. From the plots 
in (a) we can see that, as a — > 0, we recover the isotropic solution (17), since in this limit £? 4 — > and 
T 4 -> -97.4 ~ -7T 4 . 



will obtain their expressions analytically for large and small T/a. They can also be computed 
for arbitrary values of T/a by fitting the asymptotic fall-off of the numerical solution, and 
the result is shown in Fig. 3. These coefficients enter the relation between the FG coordinate 
v and the u coordinate of (8)-(9), which takes the form 

B- ■ 7./"; , 1091a 4 <r , \ - - 



v _ u--u> + { _ + _-_l„ g „J„» + 0M , (30) 

a 2 3 /£ 4 + 7J" 4 419a 4 a 4 , \ 5 7x 
U = U+ 12 V + 32256 + 16 I0g ^ +0(V ) ' (31) 

In fact, T4 and B4 are the coefficients of the 0(u 4 ) terms in the near-boundary expansion 
of the functions T and B in terms of the u-coordinate. The change to the FG coordinate v 
introduces the additional a 4 terms in the 0(v 4 ) coefficients in eqn. (28). 

Note that the logarithmic terms in (29) are independent of T4 and B^ and are given by 
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Observe too that g^ 1 is traceless, i.e. Trg( 4 ) = gV^g^ = 0. 

Again using the results from [100], the expectation value of the stress tensor is found to 

{T ij ) = ^[ 9 f-{^)gP + C fgf + ...], (33) 

1 J I < U \IV) U V U m/jM'UU 'Ml. I ' 

the solution (29) this becomes 



where the dots stand for terms that do not depend on g or x ( '■ When evaluated on 



(Tij) = diag(T, P xy , P xy , P z ) , (34) 
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with 



TV 2 / 3_ 23 2777 4 c sch , 

£ = — ^ — .F 4 Ba H a H a 

2tt 2 V 4 28 16128 96 



r> -V, 2 / 1 5 , 611 , .„.„ 

Fxy ~ 2^ V 4 >4 " 28^ 4 + 16128° " 96 a 
N 2 / 1_ 13 „ 2227 



sch 4 



, , Ta - —Ba + -^^a 4 + —a 4 . (35) 
2vr 2 V 4 28 16128 32 j v ; 

We use the notation P xy = P x = P y as a reminder that the pressures in the x- and y- 
directions are equal. Similar expressions can be obtained [100] for the expectation values of 
the operators dual to cj>, ~ Tr.F 2 , and to \-> @x ~ TrFF, from which one finds that 
(0^) 7^ and (O x ) = 0. The near-boundary expansion determines the divergence and the 
trace of the stress tensor, which are independent of the J-4 and Ba coefficients. Specifically, 
the divergence obeys 

d l (T %] ) + (O^drf® + (O x )d jX i0) = . (36) 

In our case this immediately implies the conservation of the stress tensor, d l {Tij) = 0, by 
virtue of the facts that dj(j)^ = and (O x ) = 0. The former is simply due to the fact that 
there is no source for in our solution, namely the fact that the non-normalizable mode of 
the dilaton vanishes, = 0. The latter is due to the fact that, although there is a source 
for the instanton density O x because the non-normalizable mode of the axion is non-zero, 
7^ 0, this does not induce an expectation value for O x . We thus confirm that translation 
invariance is preserved, as anticipated above. Similarly, the trace obeys 

(Ti)=A(g ( i ?,^\x i0) ), (37) 

where 

A{gf^\x^) = ^ (^diX^X^Y = Hm V=M(W>X) = ^ (38) 

is the conformal anomaly evaluated on our specific solution. This type of 'matter conformal 
anomaly' is related to short-distance singularities in the two-point function of O x and was 
first studied in [108], where it was shown that it does not renormalize in the sense that its 
coefficient is fixed relative to the normalization of the two-point function. On the gravity side 
the anomaly was first computed in [109] and is directly related to the fact that some of the 
coefficients of the logarithmic terms in the expansion above, concretely gj 4 and (f>^ , do not 
vanish. We will see that the anomaly has important consequences for the thermodynamics 
of our system. 

For our purposes it is crucial to understand how the terms in the near-boundary expansion 
transform under a rescaling of the coordinates of the form 

Xi = k x\ , v = kv' , (39) 

where A; is a positive real number. In the gauge theory this is equivalent to the rescaling of 
a and T 

a' = ka, T' = kT, (40) 
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as is easily verified from their definitions in terms of parameters on the gravity side: 

x = aZ =(ka)z', (A;T ) = _*L = i_. (41 ) 

Following [104], we note that the rescaling above leaves the FG form of the metric (22) 
invariant and transforms all the expansion coefficients homogeneously, 

9{o)ij = 9(o)ij , g[ 2 )ij = k2 9(2)ij , 9(4)ij = k 9(4)ij , ( 42 ) 
except for guuj, which acquires an inhomogeneous piece due to the logarithmic term in (23): 

0(4)y = k 4 9(4)ij + 2 k 4 log k g^ i:j . (43) 

An analogous statement holds for the transformation properties of the expansion coefficients 
of 4> an d X- It follows that the stress tensor (33) transforms as 

(Tij) = k 4 (T i:j ) + 4^ 4 log k ~g m , (44) 

namely 

N 2 n 4 

(Tijika, kT)) = k 4 (Ti^T)) + k 4 log k ^J^y ^ > ( 45 ) 

where we have made use of (32). This immediately implies that the stress tensor must take 
the form 

(T^a, T)> = a 4 Uj (£) + a 4 log ( -) ^ , (46) 



where [i is some arbitrary reference scale, a remnant of the renormalization process much 
like the subtraction point in QCD. The first and second terms transform homogeneously and 
inhomogeneously under the rescaling (40), respectively. Needless to say, one could rewrite 
the first term in a variety of forms, e.g. T 4 tij(a/T), etc. Also, one could replace log(a//i) by 
log(a/T) + log(T///), thus redefining 

^■^% + 4^2 1 °g(|) /l 'i- ( 47 ) 

In any case, the key conclusion is that the physics does not only depend on the ratio a/T, 
but on the two independent dimensionless ratios that can be built from a, T and [i. 

We now realize that the freedom in the choice of scheme associated to shifts of c sch in the 
counterterm action (21) is precisely equivalent to a rescaling of the reference scale [i. This 
follows from the fact that, as we see from eqs. (33), (32) and (46), both transformations shift 
the stress tensor by an amount proportional to A x hij. The fact that one must introduce 
a reference scale /x in order to define the theory with o / is a direct consequence of 
the conformal anomaly. There is a (limited) analogy with the situation in QCD with one 
quark flavour with M q ^ 0. The classical QCD theory has only one scale, M q . However, 
upon renormalization one realizes that a scale /x must be introduced in order to define the 
quantum theory. This scale is arbitrary, and one can trade it, for example, for the scale A QCD 
at which the coupling constant becomes strong. Similarly, in the present case the classical 
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Figure 4: The energy and pressures normalized by their isotropic values (48) as functions of a/T, 
with T ~ 0.33 in (a) and T ~ 1.1 in (b). We have chosen c sch = —1 and /i = 1. 



theory (at zero temperature) has only one scale, a, but the quantum theory depends on a 
and \i. The reason this analogy is not exact is that the scale A QCD in QCD survives even 
in the limit Af q = 0, whereas in our case the theory with a = is truly conformal: in this 
limit there is no anomaly and the scale \x disappears. Note also that, if we were interested 
in working with the theory at a fixed value of a, we could choose [i = a thereby seemingly 
eliminating the //-dependence all together. However, we will be interested in examining how 
the physics depends on a, in which case the choice above is not appropriate. For example, in 
taking derivatives of the free energy with respect to a to obtain the chemical potential, as in 
(66), a and \x must be treated as independent quantities. 

Below we will obtain analytical expressions for the energy and the pressures (35) in the 
limits of large and small T/a. For arbitrary T/a their values can be obtained numerically. 
The results for some specific cases are shown in Figs. 4 and 5. In Fig. 4 we fix T ~ 0.33 and 
T ~ 1.1 and normalize the results by the isotropic N = 4 SYM values 

o 2 N 2 T 4 7r 2 /V 2 T 4 

e°(t) = XX • = ■ ^ 

We see from Fig. 4 that in the limit a — > the energy and the pressures approach their 
isotropic values, as expected. In Fig. 5 we fix a ~ 0.34 and a ~ 2.86 and normalize the 
energy and the pressures by N^a 4 . Note that the energy and/or the pressures can become 
negative at low temperatures depending on the value of a/fj,, as is in fact apparent from (46). 
We will come back to the interpretation of these results below. 
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(a) (b) 

Figure 5: The energy and pressures divided by N^a 4 as functions of T/a, with a ~ 0.34 in (a) and 
a ~ 2.86 in (b). We have chosen c sch = —1 and = 1. 



4 Thermodynamics 
4.1 Number of D7-branes 

As anticipated in the Introduction, the solution at hand can be thought of as describing a 
uniform D7-brane charge density per unit length in the z-direction [73] . To obtain the correct 
normalization, a reminder about conventions is necessary. The ten-dimensional supergravity 
action in the string frame takes the form 

S s ^ = ^ f^\ e ~ 2VR -J2^ P n + --- ) , (49) 



2re? J v * I 2 • n! 

where 2kj = (27r) and F n denote the Ramond-Ramond (RR) field strengths of different 
degrees. In particular, F\ = dx- (The reason for the tildes will be clear shortly.) The 
dilaton ip appearing in this equation is not normalized to zero at the boundary, but in fact 
determines the string and the SYM coupling constants as g^/An = g s = e^^y. Similarly, 
the theta-angle is related to the boundary value of x as # = 27rx bdry . Given a collection of 
overlapping D7-branes (see Fig. 6 (left)), their total number is measured by the circulation of 
the axion around the branes, 

w — sk/*- (50) 

as appropriate since D7-branes are magnetically charged with respect to the axion. In this 
equation, T^ 1 = (2-7r) 7 £g is the inverse D7-brane tension, and we note that 2h 2 T U7 = 1. By 
symmetry, in this case we have x = c^N D7 /2n, with a the polar angle in the plane transverse 
to the branes. Suppose now that we distribute the D7-branes homogeneously along the z- 
direction, as in Fig. 6(right), then a similar symmetry argument shows that in this case 
X = n D7 z, where n D7 = N D7 /L Z is the density of D7-branes along the z-direction. 5 It follows 



5 Unless otherwise stated, in this paper we take Not, L z — > oo while maintaining nm finite, i.e. the 
2-direction is non-compact. One could consider instead the ^-coordinate to be periodically identified with 
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Figure 6: Distributions of D7-branes. 



that the 0-angle is given by 

6 = 2irn or z . (51) 

In this paper we found it convenient to work with a dilaton <p normalized in such a way 
that 4> hdry = 0. This can be achieved by defining 

= g s e v , F n = g s F n . (52) 

Under this rescaling the supergravity action rescales homogeneously with a factor of g~ 2 and 
takes the form 

Fn + --- I , (->>) 



where now 2n\ Q = (2tt) (-s9s- The fields appearing in this action are the ones we have used 
throughout this paper. Similarly, the D7-brane action transforms as 



-T D7 J e~^^~g + • ■ ■ 

-T D7 / e-^^-g + ■■■ , (54) 



and therefore T^ 1 = (2vr) 7 £fc/ s . We see that the coefficient a in the rescaled axion % = az is 
given by 

a = g anm = -— , (55) 

where A = = 4ng s N c is the gauge theory 't Hooft coupling. Interestingly, the factor on 

the right hand side is analogous to the one found [110, 111, 112] 6 to control the backreaction 



period L z < 00. This was the choice adopted in [91]. The fact that the action (1) (with anti-periodic 
boundary conditions in the z-direction for the fermions) flows in the IR to a 3-dimensional pure Chern-Simons 
theory was used in that reference to holographically model the fractional quantum Hall effect. 
6 Note that the convention used in [110, 111] was g Y M = 2ng a . 
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of flavour D7-branes on the geometry dual to J\f = 4 SYM plasma with quarks [96]. Note 
that in this paper we keep a fixed as N c — > oo in order to obtain a finite limit, which implies 
that the 0-angle and the D7-brane density scale as 



e 



(56) 



~ n m 



A ' 



The power of N c is expected from the fact that 9 /N c should be kept fixed in order to achieve 
a finite limit as -/V c — > oo, and also from the fact that the number of D7-branes must be of 
the same order as the number of D3-branes since they produce an 0(l)-modification of the 
geometry. The power of the 't Hooft coupling is the same as that encountered in [110, 111, 112] 
and is a strong-coupling result. 

4.2 Chemical potential 

In the 'magnetic description' in terms of the axion one cannot introduce a chemical potential 
conjugate to the number of D7-branes [113]. In other words, in this description one 
is restricted to work with a constant total number of D7-branes. In order to explore the 
complete thermodynamics of the system, we must first clarify how a non-constant number of 
D7-branes can be described by working with the field with respect to which the D7-branes 
are electrically charged. In ten-dimensional type IIB string theory this is the RR 8-form 
defined (roughly) through dC$ ~ *dx, where ★ denotes the Hodge dual in ten dimensions. It 
is easy to see from this relation that C% may be taken to have only one non-zero component 
C txy gi ...05(ii), where 9 a are coordinates on the S" 5 . This confirms our intuition that the D7- 
branes are wrapped on the S 5 , extend along the ixy-directions and are smeared along the 
z-direction. Thus from the effective five-dimensonal viewpoint the D7-branes behave as a 
2-brane charge density oriented in the xy-directions and distributed along the z-direction. 

Upon reduction on the S 5 , Cg gives rise to a 3-form that we will denote simply as C. Its 
precise definition in the effective five-dimensional theory follows from the requirement that 
its Bianchi identity be precisely the equation of motion for the axion resulting from the action 
(7) (and vice versa). Thus 



where in this and subsequent equations * denotes the five-dimensional Hodge dual with re- 
spect to the metric (8). Changing to the FG coordinate via (31), writing C = Ct xy (v) dt A dx A dy, 
and using the near-boundary expansions (28), we see that asymptotically 



The integration constant C is fixed by the regularity condition that Ctxy vanish at the hori- 
zon [114]. The most straightforward way to understand this requirement is to continue to 
Euclidean signature. In that case dt E is ill-defined at u = u H (the 'tip of the cigar') and 
hence Ct EXy must vanish there. Equivalently, in Lorentzian signature dt is ill-defined at the 
bifurcation surface of the Killing horizon, and therefore Ctxy must vanish at the horizon. 



dC = e 2rp *d X = a 



u 3 



dt A dx Ady A du 



(57) 




(58) 
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We must now pause to appreciate an important difference between a 1-form and a 3-form 
gauge potential in an asymptotically AdS§ spacetime. The time component of a 1-form 
potential falls off as 

A t = ci + c 2 v 2 + ■■■ . (59) 

In this case, up to a normalization constant, one identifies the coefficient of the leading term, 
ci, with the chemical potential, and the coefficient of the subleading term, C2, with the charge 
density. This identification is consistent with the fact that in this case *dA near the boundary 
of AdS§, whose integral measures the total charge, is proportional to C2 and independent of 
c\. In contrast, in the case of a 3-form, one must identify the coefficient of the leading 1/v 2 
term (times —2) with the charge density, as is clear from the fact that this is proportional to 
*dC ~ a and independent of C. Correspondingly, one must identify the subleading constant 
term C with the chemical potential, up to a normalization constant. As we will see below, 
the normalization 

C N 2 C 

is the correct one, in the sense that the potential defined in this way agrees exactly for all 
values of a and T with the thermodynamic quantity defined in (66), provided we choose 
c SC h = — 1- Details can be found in the next section and in Appendix D. Suffice it to note 
here that the fact that one must choose a particular value of c sch is not at all surprising. 
Indeed, a renormalized 'electric' action in which the axion had been traded in favour of the 
3-form C would still suffer from scheme-dependence associated to the possible addition of 
finite counterterms. By writing eqn. (60) we have implicitly made a particular choice of 
scheme in this electric formulation, as is clear from the fact that the rescaling (39) would 
shift C due to the logarithmic term in (58). Thus it is not surprising that the definitions (60) 
and (66) coincide only provided a particular choice of scheme is also made for the free energy 
F that enters the latter equation. We will verify below that, once this choice is made, any 
subsequent change of scheme respects this agreement. 

4.3 Fundamental thermodynamic identities 

As usual, the on-shell classical action (20) provides a saddle-point approximation to the 
corresponding thermodynamic potential through the identification 

S on _ shcll = p J d 3 xF, (61) 

where /3 = 1/T. In this case, since we are working at constant temperature and charge density, 
F(T,a) is the free energy density in the canonical ensemble (see e.g. [113]). Therefore F is 
related to the energy density through 

F = E — Ts , (62) 

and it obeys 

dF = -sdT + <S>da , (63) 
which together with (62) is of course equivalent to the first law of thermodynamics 

dE = Tds + §da . (64) 
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Note that eqn. (63) defines the chemical potential as the intensive variable conjugate to the 
charge density a, namely we have 



dF 
df 



(65) 



We will verify below that this definition of $ agrees with the definition (60). Note that a has 
dimensions of energy and that <I> has dimensions of energy/length 2 . This is appropriate since 
a measures the number of D7-branes per unit length in the z-direction, and <I> measures the 
energy cost per unit xy-area of introducing a D7-brane that extends in those directions - 
see Appendix C. 

The thermodynamic potential in the grand canonical ensemble is simply given by 

G = E-Ts-<S>a, (67) 

and it obeys 

dG = -sdT - ad<f> . (68) 

Presumably, G would equal the on-shell action computed in terms of the 3-form C dual to 
X [H3], once an appropriate renormalization procedure had been applied in such an 'electric 
formulation'. 

Extensivity of the energy E(s, a) as a function of its natural variables implies that the 
thermodynamic potentials are related to the pressures through (cf. [72]) 

F = -P xy , (69) 
G = -P z , (70) 

which implies that the pressure difference is sustained by the presence of a non-zero chemical 
potential, namely 

P z - P xy = $a . (71) 

The relation G = —P is familiar in the isotropic case. A derivation of the relations (69)-(70) 
in the anisotropic case is given in Appendix C. 

An important consistency check is the fact that all the thermodynamic identities above 
are scheme- independent, meaning that they are invariant under rescalings of fi. Consider for 
example the relation (62). From eqns. (20) and (21) it is clear that under the rescaling (39) 
F transforms as 

7V 2 o 4 

F(ka, kT) = k 4 F(a, T) + A; 4 log k -t%-~ , (72) 



48vr 2 



which immediately implies that F must take the form 



F(a,T) = a*f(?-) + a*lo g («) -g, . (73) 



Therefore 

dF dE 



TJ b \fiJ 48vr 2 

(74) 



d/j, d/j, 
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consistently with the fact that s is scheme-independent, 

as expected since it is evaluated at the horizon. In contrast, we emphasize that the chemical 
potential is scheme-dependent. This is suggested by the fact that it cannot be obtained with 
horizon information alone (since in fact regularity implies that it must vanish there) but 
must be obtained from the fall-off (58). The scheme-dependence of is also implied by its 
thermodynamic definition (66), which together with (73) implies that 

To see that this is consistent with the definition (60), we note from eqn. (58) that under the 
rescaling (39) the constant term transforms as 

C' = k 3 C - -A; 3 log A;, (77) 

where the k 3 factor comes from the transformation of dtAdxAdy. Consequently, the potential 
transforms as 

N 2 

$(ka,kT) = k 3 <S>(a,T) + Ak 3 logk — ±- , (78) 

487H 

which implies (76). In summary, eqs. (46), (73) and (76) imply that 

dE = _8Pxy = idP 1 = ad$ 
dfi dpi 3 d[i 4 d[i ' 

which together with (75) guarantee that the thermodynamic identities (62), (67), (69) and 
(70) are scheme- independent. 

Incidentally, we note that the relations above, together with the anomaly equation 

(T})=-E + 2P xy + P z = A, (80) 

imply that 

E--Ts-- A <S>a=-A. (81) 
4 4 4 

As discussed e.g. in [88, 89], 7 if the only two scales were truly the temperature and the 
chemical potential then the first law would imply that the combination on the left-hand side 
would vanish identically. The fact that this is not the case here is another manifestation of 
the presence of an additional scale (the reference scale fj,) , as implied by the anomaly. 

To close this section, we recall that the necessary and sufficient conditions for local ther- 
modynamic stability in the canonical ensemble can be written as [115] 

(dE\ rT ,(dS\ , s 

(d§\ (d 2 F\ . . 

* s U) T =UO T >o - <83) 



7 The factor of 1/4 in front of $a differs from the 3/4 factor in those references (in four spacetime dimensions) 
because here a has dimensions of mass instead of mass 3 . 
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The first condition is just the statement that the specific heat at constant charge is posi- 
tive. The second one is the statement that the system is stable against infinitesimal charge 
fluctuations. 



5 Phase diagram 



Above we reviewed and clarified certain aspects which we will now apply to the thermody- 
namics of our specific system. In order to do so we will examine the limits of high and low 
temperature. This information alone will suffice to construct the qualitative phase diagram 
shown in Fig. 9. 

In the high-temperature limit, T 3> a, [a, one can obtain approximate analytical expressions 
for the various quantities of interest by perturbing around the isotropic solution (17). The 
details can be found in Appendix D. Here we simply collect the results for the energy density, 
the pressures and the entropy density. Omitting terms of 0(a 6 ) these are 
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(84) 
(85) 
(86) 
(87) 



where E°(T),P°(T) and s°(T) are the corresponding results for the undeformed (i.e. a = 0) 
M = 4 SYM theory given in (18) and (48). The sign of the leading, 0(a 2 ) correction to 
these results is consistent with the positive slope of the E and P xy curves and the negative 
slope of the P z curve in Fig. 4 at small a/T. The chemical potential can be obtained through 
eqn. (66) and it reads 



N 2 T 2 a N 2 a 3 
$ = — h 
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(88) 



In the opposite limit, T <C a, /i, one can also obtain approximate analytical expressions for 
the various quantities of interest. The details can be found in Appendix E. Here we simply 
collect the results for the energy density, the pressures and the entropy density: 



E 



P- 



2„4 



N 2 a 
5376vr 2 
N 2 a 4 



5376vr 2 
N 2 a A 



(28c 
(28c sch 



-l + 192c int ) 
1 + 192c i; 



N 2 a A , / a 
+ 4l^ l0g 



N 2 a A , / a 



11 

3c, 



+ 



N 2 a 4 

53767r 2 ( 84c - + 109 + 576c -) + 3 ^ bg 
c eat N?a 1/3 T s/3 + • • • , 



11 

2c 



N y/3 T u/3 + . . 

- iV 2 a l/3 r ll/3 + 



(89) 
(90) 



H-^iVV/3 T ll/3 + ... )(91) 



(92) 



21 



where we have omitted terms of higher order in an expansion in T/a. In these expressions c int 
is an integration constant that could be determined numerically, and c cnt > is the constant 
defined in (19). The specific values of these constants will have no effect on the qualitative 
physics. The chemical potential is 



Note that at T = one has 8 



iV 2 a 4 

i^-P^, SP^ + P^^. (94) 



These results are in agreement with Fig. 5, where in particular we see that E and P xy have 
opposite signs in the limit T/a — > 0. The signs of E,P xy and P^ at low T/a depend on the 
value of the ratio a/\x. For example, the energy becomes negative as T — > provided a is 
sufficiently smaller than \i. We emphasize that this is not an inconsistency: It means that 
the energy can be lower than that of the undeformed M = 4 theory at zero temperature, but 
the energy is still bounded from below for all values of a, and achieves this absolute minimum 
at a scale set by fi. We will come back shortly to the meaning of the signs of the pressures. 

The coefficients of the log(a//i) terms above, both at high and at low temperature, are in 
perfect agreement with those in eqs. (46) and (76), as is the trace of the stress tensor 

N 2 a 4 

{Tl) = -E + 2P xy + P z = ^ (95) 

with eqs. (37)-(38). Also, with the explicit expressions above in hand, it is straightforward 
to verify that all the thermodynamic identities of Section 4.3 are identically satisfied in both 
limits. 

We are now ready to understand the phase diagram qualitatively. We begin by examining 
the system at zero temperature and finite density. In this case we see from eqs. (90), (91) 
and (93) that 

N 2 a\ /a\ . s 

log - , (96) 
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where we recall that < 3? / = d^/da and F = —P xy . Simple inspection of these equations 
shows that F, <J>, P z and are all positive for a S> /J>, but that as a gradually decreases 



8 In this limit the horizon disappears and the string-frame metric develops a naked curvature singularity 
deep in the IR (see Section 6 and Appendix B for details). To avoid this one can simply imagine that T is 
made very small but not zero. 
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they become negative sequentially in the order in which they are listed above, at values 
a F > a s > a Pi > a S '. Note that the precise values of 04 depend on c sc h and c- int , but that 
their ordering is independent of these constants, and that the scale of all the Gtj is set by 
/i. It follows that at zero temperature and sufficiently low densities, a < a&, the condition 
(83) is violated and therefore the system is unstable against infinitesimal charge fluctuations, 
i.e. small clumps of charge with density slightly higher than the average will grow instead of 
relaxing back to the average. At densities a 4> ' < a < a Pr the system becomes stable against 
infinitesimal charge fluctuations, and therefore we will refer to this phase as metastable. Yet, 
the system is still unstable against finite charge fluctuations, as signaled by the fact that if 
a < a Pr then the pressure of the isotropic phase is higher than that of the anisotropic phase: 9 

P z (a,T) <P°(T). (100) 

As a consequence, bubbles of the isotropic phase can form and grow inside the anisotropic 
phase, forcing a redistribution of the total charge into a smaller volume. 10 In other words, even 
if one could somehow prepare a homogenous phase of density a < a Pz , this would fall apart 
into a mixed phase consisting of 'droplets' or 'filaments' (or, more precisely, thick membranes) 
of anisotropic phase with density a = a Pz and P z = surrounded by vacuum regions with 
a = P° = 0, as depicted in Fig. 7(left). Note that the chemical potential satisfies $(a P J < 0. 

111 11 inn 

Figure 7: (Left) Inhomogeneous mixed phase in which isotropic and anisotropic regions coexist. 
(Right) Homogeneous anisotropic phase. 

This is an important consistency condition, for it means that it is indeed advantageous for 
the charge to stay together in the anisotropic region; we will comment further on this below. 
At densities a > a P% the preferred phase is a homogeneous anisotropic phase, as shown in 
Fig. 7(right), which is stable. This physics is similar to that of zero-temperature QCD at 
finite baryon density [116, 117, 118], which suggests that the phase transition from the mixed 
phase to the homogenous anisotropic phase may occur by a percolation mechanism as the 
anisotropic regions merge, again in analogy to models of restoration of chiral symmetry in 
QCD [119, 120]. We will elaborate on this similarities in Section 6. 

The dynamics of the isotropic bubbles may be complicated and will depend, among other 
things, on their surface tension. We will come back to this point in Section 6. Here we 
just note that the 'mechanical condition' (100) is exactly equivalent to the condition that a 
compression of the charge is thermodynamically preferred. Indeed, let L z be the length of 

9 In this and in subsequent equations we will not set T — explicitly because, as we will see, this instability 
extends to T > 0. In order to specify to T — one must simply remember that P°(T) — 0. 

10 Note that no charge redistribution can occur in the xj/-directions, since the branes are extended along 
these directions. 
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the z-direction, which can be taken to infinity at the end. We wish to fix the total D7- 
brane charge a (we follow the notation of Appendix C) and compare the free energies of the 
following two configurations (see Fig. 8): (i) a configuration in which the charge is distributed 



(l-e)L s eL z 



Figure 8: Redistribution of D7-brane charge along the z-direction. 



homogeneously over the entire length L z with density a = d/L z , and (ii) a configuration in 
which a small fraction eL z of the total length is occupied by the isotropic phase and the rest 
is occupied by the anisotropic phase with the total charge distributed homogeneously but 
with a higher density 

a' = a ~ a + ae. (101) 
(.1 — £)L Z 

The free energy (per unit xy-area) of each configuration is 
F {1) = L z F(a, T) 

F m = L z (l-e)F(a + ae,T) + eL z F°(T) 



'(«) 

~ L z F(a, T) + eL z F°(T) - eL z ( F(a, T) - dF ^l T l ' a 

\ da 

= L z F(a, T) - eL z P°(T) + eL z [p xy {a, T) + $(a, t)aj 

= L z F(a,T)-eL z P°(T) + eL z P z (a,T), (102) 

where we have used eqs. (66), (69) and (71). We therefore conclude that Fu^ < Fr^ if and 
only if (100) is satisfied. 

In order to understand how the physics above extends to or is modified at T > 0, let us 
examine the limit T 3> a. Before addressing the behaviour under charge fluctuations, we 
briefly establish that there are no thermal instabilities. Indeed, ignoring terms of 0(a 4 ) in 
(84) we see that the specific heat (82) is positive at high temperatures. The first temperature- 
dependent correction at low temperatures in eqn. (89) also yields c a > 0. We see from the 
positive slope of the red, continuous curve in Fig. 5 that this extends to all values of T/a. 
In fact, not only the energy but also the pressures are monotonically increasing functions of 



'Note that the N c — > oo limit in which we are working ensures that we are in the thermodynamic limit 



even if the total volume V = L x L y L z is finite. 
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Figure 9: Qualitative phase diagram. The orange, dashed-dotted curve is defined by the condition 
F(a,T) — F°(a,T). The magenta, dotted curve is defined by the condition $(a,T) = 0. The blue, 
continuous curve is defined by the condition P 2 (a,T) = P (T). The dashed, red curve is defined by 
the condition Q'(a, T) = 0. In Zones I, II and III the system is in a homogeneous phase, whereas it is 
in a mixed phase in Zones IV and V. See the main text for more details. 



T for fixed a. This implies that the speed of sound (squared) in the xy- and z-directions is 
positive: 

4=(^) # >o, «2=(^).>o. ao3> 

Returning to charge fluctuations, we see from (88) that 

/<9<I>\ N 2 T 2 

at sufficiently high temperatures, demonstrating that the local stability condition (83) is 
violated in this limit. We also note from eqn. (86) that the leading 0(a 2 ) correction to P z is 
negative, meaning that we are again in the situation described by eqn. (100). 

With this information we can now draw the qualitative phase diagram — see Fig. 9. There 
are four curves denned (from right to left) by the conditions F(a, T) — F°(a, T) = (orange, 



25 



dashed-dotted curve), $(a, T) = (magenta, dotted curve), P^a, T) - P°(T) = (blue, solid 
curve) and &(a,T) = (red, dashed curve). These four curves intercept the horizontal, zero- 
temperature axis at the points a» defined above. The slopes at T = of the first two curves 
are positive, whereas the slopes of the other two are negative. This can be seen by examining 
the leading temperature correction to the zero-temperature results (98)-(99). Consider for 
example the chemical potential. Let (nT,n a ) be the tangent vector to the curve $(T, a) = 
at the point (T,a) = (0, o*). By definition the derivative of $ along the curve vanishes, 
namely 

d<S> d<S> 







nT df +Ua d^ 



(105) 



Obviously at (T, a) = (0, a $ ) we have tit > 0, and we see from eqn. (93) that d^/dT < 0. 
Also, we know that dQ/da > since $ goes from negative to positive across the curve in the 
direction of increasing a. It follows that we must have n a > 0. A similar argument applies 
to the other three curves. 

Note that the four curves must 'extend to infinity' in the sense that they must contain 
points for which both T and a are arbitrarily large. To see this, consider again the chemical 
potential. This is negative in the limit T ^> a, \x but positive in the limit a>T,/i, Thus if we 
start at any point (a, T) at which is negative and we increase a keeping T fixed, eventually 
we reach the regime in which a>T,/j where must be positive. Thus a horizontal, constant- 
T line always crosses the curve $ = no matter how large T is. Similarly, if we start at 
a point (a, T) where $ is positive and we increase T keeping a fixed, eventually we reach 
the regime in which T ^> a, fi where <!> must be negative. Thus a vertical, constant-a line 
always crosses the curve $ = no matter how large a is (specifically, provided it is larger 
than a s ). An analogous argument applies to each of the other three curves. Of course, from 
the asymptotic information at large and small a/T we cannot deduce the detailed form of 
the curves at arbitrary values of a/T; for example, we cannot exclude that some curves cross 
each other. However, we have verified the qualitative form of the phase diagram shown in 
Fig. 9 with our numerical solution. In particular, we have verified that the curves do not 
cross. 

Returning to the phase diagram, we conclude that there are five distinct zones: 



I 

II 
III 

IV 

V 



F(a,T) > F°(T) 
F(a, T) < F°(T) 
F(a, T) < F°(T) 
F(a, T) < F°(T) 
F(a, T) < F°(T) 



$(a,T)>0, P z (a,T) > P°(T) 

$(a,T)>0, P z (a,T)> P°(T) 

$(a,T)<0, P z (a,T)> P°(T) 

$(a,T)<0, P z (a,T) <P°(T) 

$(a,T)<0, P z (a,T) <P°(T) 



&(a,T)>0, 
$'(a,T) >0, 
&(a,T)>0, 
$'(a,T) >0, 
&(a,T) < 0. 



(106) 



In Zones I, II and III the thermodynamically preferred configuration is a stable homogeneous 
phase. In contrast, the homogeneous phase is metastable in Zone IV and unstable in Zone V. 
In these zones the preferred configuration is a mixed, inhomogeneous phase in which multiple 
regions of the isotropic and the anisotropic phases coexist with each other, as in Fig. 7. The 
local charge density in each of the anisotropic regions is a ani > a and the pressure exactly 
equals that of the isotropic phase at the same temperature, P z (a a _ ni ,T) = P°(T). The pair 
(a ani ,T) therefore lies on the blue, continuous curve of the phase diagram. Note that the 
chemical potential need not be, and in fact is not, the same in the isotropic and anisotropic 
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regions: $ vanishes in the isotropic phase and is negative in the anisotropic phase for points 
that lie on the blue curve. As mentioned above, this is an important consistency condition for 
the coexistence of the two phases. In a system in contact with an infinite charge reservoir the 
chemical potential will eventually equilibrate to the same value everywhere. However, in the 
system under consideration, once an anisotropic region is surrounded by isotropic charge-less 
regions, the anisotropic region cannot draw additional charge from the latter. In this case 
consistency only requires that the chemical potentials verify <I> ani < $ iso , since otherwise the 
charge would escape from the anisotropic to the isotropic regions. 

If we vary the temperature of a mixed-phase configuration adiabatically, P° (T) will change 
smoothly and the volume occupied by the anisotropic regions will adjust so that the pair 
(a ani ,T) stays on the blue curve. For example, if we start with a density and temperature 
(a,T) corresponding to the point '1' in the phase diagram, the corresponding pair (a ani ,T) 
describing the anisotropic regions will be point '2'. As we lower the temperature, the pair 
(a, T) moves down along a vertical line, whereas (o ani ,T) moves down and leftwards along 
the blue curve. The two trajectories meet at point '3'. At this point the two densities agree, 
a = a an i) and thus the anisotropic phase fills out the entire space, i.e. we have reached the 
homogeneous phase. As we keep lowering the temperature we stay in the homogeneous phase. 
Because of the form of the blue curve, it is clear that if we start with a density a to the left 
of the turning point of the blue curve then as we vary the temperature we never reach the 
homogeneous phase. Similarly, for an intermediate range of initial densities between the 
turning point and a Pz , we may enter and subsequently leave the homogeneous phase. 

6 Discussion 

In this paper we have studied the finite-temperature generalization of the type IIB super- 
gravity solution of [73] dual to an anisotropic deformation of N = 4 SYM. The solution 
obeys AdS^ x S 5 boundary conditions, as appropriate for a marginally relevant deformation, 
possesses a regular anisotropic horizon, and is completely smooth on and outside the horizon. 
These features guarantee that it is solidly embedded in string theory. The deformation of the 
gauge theory is induced by a position-dependent ^-parameter (see eqn. (1)), which in the dual 
gravity description can be understood as resulting from a density of dissolved D7-branes. 

The fact that the D7-branes are completely wrapped on the S* 5 leads to a great technical 
simplification, because it implies that the 50(6) symmetry of the undeformed theory is 
preserved. This means in turn that none of the Kaluza-Klein (KK) modes on the 5 5 are 
excited, and therefore that the solution can be found by working directly with five-dimensional 
supergravity with a negative cosmological constant coupled to a few matter fields — just the 
axion and the dilaton. 

At zero temperature the IR limit of the solution is given by eqs. (144)-(145). As explained 
in Appendix B the string-frame metric possesses a naked singularity deep in the IR, i.e. at 
r — > in the coordinates of (144) or at u — > oo in terms of the coordinate used above. The 
metric in the Einstein frame is also pathological. As pointed out in [121] and elaborated 
upon in [122], although all curvature invariants of the Einstein- frame metric are finite, tidal 
forces diverge as one approaches the r = hypersurface. In contrast, at finite temperature 
the problematic region is hidden behind a horizon. It is remarkable that in this case the 
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solution takes into account the full back-reaction of the D7-branes but exhibits no pathologies 
regardless of their number (density), which scales as n D7 ~ N c /X (see eqn. (56)) and can be 
larger than 24. This provides a completely explicit counterexample (at least in the presence of 
a negative cosmological constant) to the somewhat extended belief that 24 is an upper bound 
on the possible number of D7-branes. In particular, the dilaton is finite everywhere. In this 
respect it is interesting to note that in the IR Lifshitz geometry (144)-(145) the dilaton grows 
without bound towards the UV (r — > oo), but this growth is tamed by UV-completing the 
Lifshitz geometry with an asymptotic AdS region. We note that the parameter controlling 
the backreaction of the D7-branes is not just the ratio of their number (density) to the 
number of D3-branes, but this ratio enhanced by a power of the 't Hooft coupling, An D7 /iV c , 
as in the case [110, 111, 112] of flavour D7-branes dual to quark degrees of freedom in the 
N = 4 SYM plasma . In the latter case, however, inclusion of the D7-branes' backreaction 
does lead to a singularity in the UV part of the geometry [123, 112]. The dual gauge theory 
interpretation is that the addition of matter degrees of freedom to a conformal field theory 
makes the /3-function positive and thus leads to a Landau pole in the UV. In contrast, in the 
present case we are not adding new degrees of freedom because the D7-branes do not extend 
to the AdS boundary, so from this viewpoint it is perhaps not surprising that the gravity 
solution does not exhibit any pathologies in the UV. 

It would be interesting to generalize our construction to dissolved branes of other dimen- 
sions. One could certainly consider other branes partially wrapped on the S 5 in such a way 
that their non-compact directions coincide with some of, but not all, the gauge theory direc- 
tions, since we want them to break isotropy. However, this may be technically challenging if 
the branes break too much of the SO(6) symmetry of the sphere, since in this case the KK 
modes will become excited. In this respect it might be useful to smear the branes on the 
S 5 as to preserve the maximum amount of symmetry possible, along the lines of [124, 125] 
(see [126] for a review). An additional desirable requirement is that the branes responsible 
for the anisotropy do not change the UV properties of the theory, since this would destroy 
the AdS asymptotics, thus obscuring the microscopic nature of the dual theory (if it exists 
at all). Presumably, this requirement will exclude lower dimensional branes, as illustrated 
by the dissolved Dl-branes and strings of Refs. [127, 10], which are dual to non-commutative 
deformations of the J\f = 4 SYM theory that do change the UV physics. It would also be in- 
teresting to generalize our construction to gauge theories in different numbers of dimensions. 
A step in this direction was already given in [73] , where an anisotropic Lifshitz- like geometry 
based on the D4-D6 system was constructed. 

Our main focus in this paper was the study of the thermodynamics of the system. For 
this purpose we obtained the holographic stress tensor particularizing the results of [99, 
100] , and showed that its transformation under the rescaling (T, a) — > (kT, ka) contains 
an inhomogeneous piece (see eqn. (45)) induced by the presence of a non-zero conformal 
anomaly A = iV c 2 a 4 j 48ir 2 (see eqn. (38)). This implies that the quantum theory depends on 
an additional scale fi, a remnant of the renormalization process, analogous to the subtraction 
point in QCD. The main result of our study is the phase diagram of the system in the density- 
temperature plane, shown in Fig. 9. As noted above, in the strict T = limit the horizon 
disappears and the string-frame metric exhibits a naked curvature singularity deep in the IR. 
In this region we expect stringy a'-corrections to become important and therefore some of 
the details of the phase diagram in a thin horizontal slice that includes the T = axis may 
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Figure 10: Cartoon of the formation of a bubble of isotropic phase. 



be modified. 

We identified zones of the phase diagram in which the homogeneous phase is metastable 
or unstable. All the instabilities are related to the tendency of the D7-brane charge to clump 
together. This follows from the fact that if the charge density in a given region is too small, 
then the pressure in the z-direction in that region becomes lower than the pressure of the 
isotropic phase at the same temperature. Once this happens, 'bubbles' of the isotropic phase 
can form and grow, thus forcing a compression of the charge density into smaller regions, as 
depicted in Fig. 7. We verified that this 'mechanical conclusion' was precisely equivalent to 
the condition that the free energy of the system be minimized. Note that in our analysis we 
assumed that the total number of D7-branes was fixed, but that their local number density 
was not. In other words, we considered inhomogeneous phases, which is of course natural 
from the viewpoint of the thermodynamics of the gauge theory. In the context of AdS/CFT 
one could choose to fix the D7-brane density at each point, since this is specified by the 
non-normalizable mode of the axion. However, one might suspect that the instability relies 
on IR near-horizon physics. If this is the case, then it is possible that even if the axion 
non-normalizable mode held fixed an instability may still arise through the tendency 

of the instanton density O x ~ Tr FF to develop an expectation value. 

Analyzing the formation and dynamics of the isotropic bubbles would require consideration 
of many effects, in particular the bubble surface tension. This is certainly beyond the scope 
of this paper, so we will just make two observations. First, we note that the isotropic bubbles 
can form even if the xy-directions are non-compact. A cartoon of how a typical bubble 
may start is shown in Fig. 10. The key point is that the bubble can form with a finite size 
in all directions, thus requiring a finite fluctuation, and subsequently expand, in particular 
in the xy-directions. Second, there should exist domain-wall solutions in the gravitational 
description which interpolate along the z-direction between regions with (1^/0 and d\ = 0, 
provided the pair (T, a) in the former region lies on the blue curve of the phase diagram. These 
solutions would be dual to the inhomogeneous configurations depicted in Fig. 7, and they 
would allow for an estimation of the surface tension of an anisotropic bubble. Presumably 
the instanton density (0 X ) ~ (TrFF) would be non-zero in the interface region. 

An interesting extension on the gravity side would be to consider time-dependent solutions. 
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The focus of this paper has been the static anisotropic solution, but one could now view this 
solution as an initial condition for a time-dependent solution. Specifically, one could look 
for a gravity solution such that the non-normalizable mode of the axion satisfies d\ 7^ for 
t < to and d\ = for t > to- Presumably this would allow one, for example, to compare the 
isotropization and the thermalization times for the resulting solution at t > to. 

The physics displayed by our system in Zones IV and V is similar to that encountered 
in QCD at low temperature and finite baryon density [116, 117, 118]. In that case one 
finds that the pressure of a chirally broken homogeneous phase with density lower than a 
critical density no is negative. 12 This indicates an instability of the system towards the 
formation of 'droplets' of higher density no (interpreted as nucleons) in which P = and 
chiral symmetry is restored, surrounded by empty space with n = and P = 0. In our case, 
the role of the chirally restored phase is played by the anisotropic phase, the analogue of no 
is a Pz , and the 'droplets' correspond to the 'filaments', or more precisely thick membrane-like 
regions, of non-zero D7-brane density. These similarities suggest that the transition from 
the mixed phase to the homogeneous phase may occur via a percolation mechanism as the 
anisotropic regions merge, in analogy with some models of restoration of chiral symmetry in 
QCD [119, 120]. Note however that the objects that percolate in our case have infinite extent 
in the rry-directions. In this respect a closer analogy might exist with suggested models of 
deconfinement based on percolation of extended colour strings [128, 129]. The similarities 
between our system and QCD extend to sufficiently low temperatures, but not beyond: In our 
case the zones where the mixed phase exists (Zones IV and V) contain points with arbitrarily 
large temperatures and charge densities, whereas in QCD the mixed phase exists only in a 
low-temperature, low-density region of the phase diagram 

The instabilities that we have uncovered are reminiscent of instabilities of weakly cou- 
pled plasmas in the presence of anisotropies, and we will now comment on some superficial 
similarities and differences between both cases. However, we emphasize from the start that 
addressing to what extent, if any, the instabilities encountered here can be viewed as a strong- 
coupling counterpart of weak-coupling instabilities is beyond the scope of this paper, and we 
make no claim to that effect in either direction. The ensuing comments should thus be taken 
merely as suggestive observations. 

Somewhat pictorially, the main similarity is the tendency to 'filamentation'. Very roughly 
speaking, in weakly coupled plasmas one can think of this as the tendency of similarly oriented 
currents to cluster together [40, 41, 42, 43, 44] (see e.g. [47] for a more precise discussion), 
as sketched in Fig. 11. We emphasize however that there are also obvious differences. In 
a weakly coupled plasma the anisotropy is 'dynamical' in the sense that it arises from the 
momentum distribution of the particles that compose the plasma, and it will disappear as 
the plasma eventually isotropizes and thermalizes. In contrast, in our case the plasma is 
static and intrinsically anisotropic because of the presence of dissolved extended objects. Yet 
the existence of instabilities in weakly coupled plasmas may be discovered by splitting the 
degrees of freedom into 'hard' and 'soft' modes [42, 43, 45, 46, 47, 48, 49, 50, 51], where the 
former carry momenta of order T and the latter of order g^mT. The hard modes are then 
treated as particles whose distribution function f(t, x, p) is often taken to be static (over an 
appropriate time scale) and homogeneous in position space but anisotropic in momentum 

12 Except in a tiny region of very small densities. 
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Figure 11: In weakly coupled plasmas parallel currents tend to cluster together. 



space, /(i,x, p) = /am S o(p)- The soft degrees of freedom are then treated as fields in this 
'hard background', and their dispersion relations reveal some unstable modes. Although in 
strongly coupled plasmas with a gravity dual there is in general no quasi-particle description 
and no separation between hard and soft degrees of freedom (see e.g. [8] and references 
therein), our gravity solution may be viewed in a similar spirit to the extent that it provides 
a homogeneous, static, anisotropic and unstable background (for certain values of (a, T)). 

It would be interesting to explore whether this analogy can be made more precise. A 
first necessary step would consist of a comparison of the degree of anisotropy in our solution 
and that of weakly coupled plasmas. In the description of the latter one often considers 
an anisotropic distribution that can be obtained by stretching or contracting an isotropic 
distribution [46, 48], i.e. one chooses 



(p) °c / iso ( VP 2 + £(P • n 



(107) 



where n is a unit vector in the direction of the anisotropy (the z-direction in our case) , 

(P 2 t) n 



£ 



(108) 



is an adjustable parameter specifying the degree of anisotropy, and (pj,) and {p\ ) are the 
average particle momenta in the directions transverse and longitudinal to the anisotropy, 
respectively. The isotropic case corresponds of course to £ = 0, and the initial stages of a 
heavy ion collision are expected to correspond to £ > [130]. In order to compare with the 
anisotropy of our solution, we note that £ is related to the anisotropy defined in terms of the 
pressure difference 

p 

A = ~p - 1 (109) 
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It is interesting to note that in our case, if one considers at least the metastable homogeneous 
phase, then A actually spans the entire range A £ (— oo, oo), part of which (A < —1) cannot 
be mapped to any value of £ > — 1. The simplest way to see this is to consider approaching 
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the blue line in the phase diagram of Fig. 9 along the T = axis. At the point (T, a Pz ) we 
have P z = and, from eqn. (94), P xy > 0. Thus as we approach the curve from the left 
A —> — oo, whereas as we approach it from the right A — > oo. We also note from eqns. (85) 
and (86) that for a small anisotropy (T 3> a) 



whereas on the magenta curve we have $ = and therefore, from eqn. (71), A = £ = 0. 

It is interesting to mention that a set of intrinsically anisotropic hydrodynamic equations, 
possibly including dissipation, can be written down for the hard modes [29, 30, 31, 32, 33, 34, 
35, 36, 37]. These equations can be derived, for example, by taking moments of the Boltzmann 
equation with an anisotropic distribution function f(t, x, p) [33, 35], and might provide an 
effective description of the system formed in the ultrarelativistic collision of two heavy ions 
in the early stage r < r iso after the collision in which the anisotropy is large. In these 
equations the magnitude of the anisotropy £ is allowed to vary over space and time, but the 
orientation of the anisotropic direction is fixed. In the absence of dissipation, the generalized 
equations in which fluctuations of the anisotropy direction are also allowed have been written 
down in [72] . In the latter case the anisotropy is induced by the presence of a p-dimensional 
charge density whose local orientation determines the anisotropic direction. Our solution may 
be viewed as the static (thermodynamic) limit of such a hydrodynamic theory with p = 2. 
Perturbing our gravity solution and applying the fluid/gravity correspondence along the lines 
of [64] it should be possible to derive the corresponding anisotropic hydrodynamic equations 
to (in principle) any desired order in the derivative expansion. Of course, the form of the 
dissipative terms at every order in the expansion can be written down on general grounds, 
but the derivation based on perturbations of our gravity solution would also predict specific 
values of the hydrodynamic coefficient of each term, as in [132, 64, 133]. Note that, unlike in 
Ref. [132], in our case the anomaly already plays a role in the thermodynamics of the system, 
and therefore it will presumably also play a role in the hydrodynamics even at the lowest 
order. 

The nature of some of the instabilities of our solution is also reminiscent of the phenomenon 
of cavitation. In an ordinary flowing liquid, cavitation is the formation of bubbles of vapour 
in regions in which the pressure of the liquid drops below its vapour pressure [134] . Cavitation 
has been proposed [135, 136, 137] as a mechanism that would lead to clustering and fragmen- 
tation of the QGP into droplets that would subsequently evaporate, thus providing a new 
scenario for how hadronization is achieved. In this context the analogue of vapour pressure 
is played by the pressure of the vacuum, P = 0, so cavitation would presumably occur if vis- 
cosity corrections to the pressure of the plasma become large enough to drive it negative. At 
that point bubbles of vacuum can form causing a fragmentation of the plasma ball. The pos- 
sibility that this takes place in the QGP has been investigated in [137, 138, 139]. The result 
depends of course on the values assumed for the transport coefficients of the QGP, in parti- 
cular on the bulk and shear viscosities. A study using holographically determined coefficients 
has appeared recently [140]. Returning to our solution, we see that the role of the vapour 
pressure is played by the pressure of the isotropic phase: the homogeneous anisotropic phase 
becomes unstable against the formation of isotropic bubbles precisely when the anisotropic 
pressure drops below the pressure of the isotropic phase at the same temperature (eqn. (100)). 
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We emphasize, however, that in the case of the mechanism proposed in [135, 136, 137] the 
pressure drop is due to a dynamical effect, namely to the viscosity corrections to the thermo- 
dynamic pressure that result from the expansion of the plasma. In contrast, in our case this 
is a static effect presumably resulting from the interaction of the D7-branes dissolved in the 
plasma. Yet, the similarity is intriguing and it would be interesting to investigate whether 
our solution can be viewed as a static toy-model for the dynamical situation. 

Before closing our discussion of instabilities, we explain why an instability studied in 
[73] does not directly apply to our setup. The authors of [73] analyzed the spectrum of 
supergravity fluctuations at zero temperature around the Lifshitz geometry (144) (in the 
Einstein frame). By KK reducing on the S" 5 , they determined the mass spectrum of the 
different modes, and by analyzing the corresponding Laplace-like equations in the metric 
(144), they conjectured a stability bound on the values of the masses, which can be thought of 
as a generalization of the usual AdS Breitenlohner-Freedman (BF) bound [141, 142, 143, 144]. 
They concluded that a few modes violate the bound, indicating an instability. The key point 
is that this analysis applies only when the solution is exactly Lifshitz everywhere, but ceases to 
apply once this solution is UV completed with an AdSs geometry. In fact, the modes identified 
in [73] as unstable with respect to the 'Lifshitz BF bound' are stable with respect to the AdS 
BF bound. Yet, it would be extremely interesting to investigate the remaining possibility 
that these modes 'want to condense' in the IR, where they are effectively unstable, in a 
way similar to that observed in the condensation of neutral scalar fields [145] in holographic 
superconductors [146, 147, 145]. Since the modes in question possess angular momentum on 
the S 5 , this would lead to a spontaneous breaking of the SO(6) symmetry of our solution at 
low temperatures. 

Throughout this paper we have implicitly assumed that the xy-directions are not compact. 
Suppose now that one of them (say x) is periodically identified with period L. If periodic 
boundary conditions are imposed on the fermions around the circle then the physics remains 
as we have described above. However, if antiperiodic boundary conditions are imposed, then 
we know [148] that the undeformed N = 4 SYM theory undergoes a thermal phase transition 
at a critical temperature T c = 1/L between a confined and a deconfined phase. In the 
gravity description this corresponds to a Hawking-Page-like phase transition [149] between 
geometries without and with a black hole, respectively. This transition extends to all values 
of a at exactly the same critical temperature because the Euclidean metrics in the two phases 
are related by a simple relabeling of the t, x coordinates. The high-temperature geometry 
is dual to the deconfined phase and is given by (8), whereas the low-temperature metric 
describes the confining phase and is given by the double analytic continuation of (8): 

ds 2 = £Zl ( _ d f2 + JT Bdx 2 + dy 2 + ndz 2 + ^\ 

u z \ F J 

Clearly, the Euclidean metrics corresponding to the Lorentzian metrics (8) and (113) are 
related by the simple exchange t o x (cf. [91]). Since the free energy density is simply 
obtained by evaluating the Euclidean action on the corresponding Euclidean solution, it 
follows that the free energy density F in the confining phase is obtained from the free energy 
density in the deconfined phase F by exchanging j3 -H- L. By symmetry, this implies that 
the phase transition must happen exactly at the point L = f3 regardless of the value of a. 
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Figure 12: Modification of the phase diagram at temperatures T < 1/L if the x- or y-dfrections are 
periodically identified with period L. 



At leading order in the large- iV c expansion F is independent of L, and therefore F and its 
derivatives <£> = dF/da and d 2 F/da 2 are temperature- independent. 14 The same is true for 
the pressure in the z-direction P z , which is obtained from P z by the exchange (3 L. Since 
the free energy, its a-derivatives and the z-pressure are continuous across the phase transition 
(though their T-derivatives are not), we conclude that the phase diagram in the (a, T)-plane 
'freezes out' below T c , meaning that the four curves in the diagram extend from T = T c to 
T = in a temperature-independent way, as shown in Fig. 12. 

We plan to report elsewhere on the effects of the anisotropy on several observables. For 
example, anisotropy effects on heavy-quark energy loss in a weakly coupled plasma have 
been studied in [150, 151, 152], and effects on momentum broadening, in particular on the jet 
quenching coefficient q, in [153]. In the isotropic case, the drag force on a heavy quark in a 
strongly coupled plasma has been considered in [38, 39, 154, 155, 156, 157, 158, 159, 160], its 
energy loss through Cherenkov-emission of heavy quarkonium mesons in [161, 162], and jet 
quenching in [163, 164, 165, 166, 167, 168, 169, 170]. Recently the drag force was also studied 
in an anisotropic background in [171] and in a Lifshitz background in [172]. Extending the 
calculations of these references to our anisotropic solution will allow a comparison with the 
results for weakly coupled anisotropic plasmas. 

Work is also in progress on the anisotropy effects on heavy quarkonium physics, which at 
weak coupling have been studied in [173, 174, 175, 176]. In our context, consideration of 

13 The masses of KK states around the x-circle depend on L, and so these states will give an //-dependent 
but iV c -subleading contribution when loop effects are considered. 

14 Incidentally, note that this immediately implies that the entropy density vanishes at leading order in the 
confined phase, § — dF/dT = 0, which is consistent with the absence of a horizon in the metric (113). 
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dynamical quarks and mesons requires the introduction of flavour branes [96] wrapping an 
iS 3 inside the S 5 . The meson spectrum and the meson dispersion relations can be calculated 
by extending the analysis of [97, 111, 177]. It is also interesting to understand the effects of 
the anisotropy on the meson dissociation phase transition uncovered in [178, 111, 179, 180, 
181]. One conclusion that can be anticipated is the fact that the in-medium meson limiting 
velocities [111] in the xy- and in the z-directions will differ, since these are respectively given 



evaluated at the tip of the flavour D7-branes at the point at which the dissociation transition 
takes place. 

A third set of observables which can be extracted from our solution is that related to 
photon and dilepton thermal production by the plasma. At weak coupling these have been 
calculated both in the absence [182, 183, 184] (see [185] for a review) and in the presence 
of anisotropy [186, 187, 188, 189]. Corresponding studies at strong coupling in the isotropic 
case include [190, 191, 192, 193], and we are currently extending these calculations to the 
anisotropic case. 
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A Derivation of the solution 

This appendix contains a detailed description of the procedure we have followed to derive 
our solution. 

The starting point is type IIB supergravity in the string frame (see for example [194]). 
Following Ref. [73] (whose relevant details we review in Appendix B), we seek a solution 
in which only the metric, the axion, the dilaton and the RR five-form are non-trivial. The 
solution can be consistently derived from the equations of motion associated to the truncated 
action 



by 




and 




(114) 



S 




(115) 



where M = 



9 and F\ = d\ is the axion field-strength. 



35 



The dilaton equation of motion resulting from this action is 

R + Ag MN (V A fVjv0 - d M <t>d N 4>) = , (116) 
and the Einstein equations are 

Rmn + 2V m Va^ + ^g M Ne^dp X d P x - \^ (f m F n + ±-F M abcdF n abcd ^ = 0. (117) 

The forms obey the following equations of motion, Bianchi identities, and self-duality con- 
straint 

d * Fx = = d * F 5 , dFi = = dF 5 , *F 5 = F 5 , (118) 

where * is the ten-dimensional Hodge dual. To model an anisotropic plasma with an asym- 
metry between, say, the xy- and the z-directions, we assume the following Ansatz for the 
string-frame metric 

1 / du 2 \ 

d s 2 = — [ -TB dt 2 + dx 2 + dy 2 + Udz 2 + — J + Z dVi% . (119) 

Without loss of generality (see Section 2), in this equation we have set to unity an overall 
constant L with dimensions of length. J 7 , B, H, and Z, as well as the dilaton 0, are taken to 
be functions of the radial coordinate u only. J7g5 is the metric of a unit five-sphere (or of any 
other Einstein five-manifold). In this parametrization the boundary is located at u = 0. As 
explained in the main text, we have used reparametrization invariance to fix g xx and g yy as 
in (119). Having done this, we cannot eliminate B in general, although we still have a scaling 
symmetry in the time coordinate that can be used to set B hdiy = 1. (Here and in the main 
text we use the subscript 'bdry' to denote the value of a field at u = 0.) Similarly, there is 
a scaling symmetry in z that we can use to set T-L hAvy = 1. J 7 is a 'blackening factor' that 
allows for the introduction of a black hole in the geometry. It vanishes at the position of the 
horizon, u = ii H - 

We take the (magnetic part of the) RR five-form to be proportional to the volume form of 
the five-sphere 

F 5 = a(tt s s +*Q S5 ) . (120) 

Equations (116)-(117) imply that at the boundary we must have a = 4e _(?ibdry . As explained 
in Section 4.1, we can choose <j) bdiy = without loss of generality, which corresponds to a = 4. 
Following [73], we choose the axion to be linear in the boundary coordinate z 

X = az, (121) 

and we simplify our Ansatz further by setting 

U = e-*, Z = e^. (122) 

Clearly, (120)-(121) satisfy (118). 

With the choices above the sphere part of the space factorizes in Einstein frame, and the 
geometry becomes a direct product of the sphere and an asymptotically locally AdS space. We 
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can therefore reduce the ten-dimensional action (115) to the five-dimensional axion-dilaton 
AdS gravity in (7), with the five-form flux giving rise to the cosmological constant. The 
metric (119) also reduces (in Einstein frame) to (8). 

Notice that the asymptotic equations of motion also imply that 

•F bdry = 1/ -^bdry • (123) 

With our choices 4> hdiy = (which imposes F^ dTy — %bdry — Zhdvy — 1) and ^bdry — 1, the 
metric (119) induces on the boundary a flat metric in which the t,x,y,z coordinates are 
canonically normalized. 



A.l Equations of motion 

We obtain now the explicit expressions for the equations of motion that we are ultimately 
going to solve. 

We start with some counting. The Einstein equations (117) are diagonal, as a consequence 
of the fact that the sphere is Einstein and that the warp factors in the ansatz depend only 
on one coordinate, u. We have one equation from (116) and five equations from (117), which 
are linearly independent (the x- and y-directions give one independent equation, and so do 
all the sphere-directions). Out of these equations, five turn out to be second-order and one 
is first-order, namely a constraint. We have only three unknown functions to solve for: (p, £>, 
and F. 

What we did in practice is the following. The six independent Einstein and dilaton equa- 
tions we started with read 

3 * , 8e"*/ 2 8F AFB' FB' 2 5F' 



e ■ 2 

~^T a ^ + u 2 uB 2B 2 



u 



3B'F' hFcp' hTB'ti 5F'<p' , FB" _„ 
+ 2B + 2u AB 4 + B + ^ ' [ ' 

e^u 2 4 4e~^ 2 B' T' 50' ,^ rS 
AF u uF 2B F 4 

^ 8 8e-^ 2 B' 2F' 110' BV F'S 5</>' 2 

-TT^a H k-=- H ^ — — = : Y<p , ( 126] 

2F u 2 u 2 F uB uF 2u 2B F 4 

e 3 * 2 16 16e-<^ 2 2B' B' 2 10 J 7 

^ a u 2 ^ + ^6 + tf 2 + 

VBF d 9d 2 2B" 2F" 

—BF- + u + ^F--^r--B--^ + 5<l) ' (12?) 

e 3 ^ 2 3^' By F'ti 5<A' 2 J/; , , 

~F a + ' (128) 

20e"^/ 2 20J" 4F8' J"fi' 2 8J 7 ' SB'J 7 ' 



u 2 u 2 



+ u£ + 2£ 2 + u 2£ 



177£ 5JV 17^0 /2 TB" hT4>" „ 
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where primes denote derivatives with respect to u. We used the last two equations, (128)- 
(129), to solve for 4>" and B" and plugged the results into the three other dynamical equations, 
(124) and (126)-(127), to obtain three constraints: 

, 6 2 12e-*/ 2 127 37' 37c// SFfffl 5F'<p' „ , /2 . , 

= 2e 3 V + = 5- + + — — + 70' 2 , 130 

u 2 u 2 u 2u AB 4 

e^u 2 8 8e-*/ 2 B' 27' hep' . , 

0= — a 2 -- + ^^ + - + — 131 
IJ- u uJ- B J- 2 

2 24 24e~^ 2 6B' 67* 18^ 50^' 57V , 2 
U " 7° + u 2 u 2 7 uB uF + « " 2£ " 7 +20 ' 

Only two of the four constraints (125) and (130)-(132) are independent. We can use for 
example (130) and (131) to solve for 7' and B'. Plugging the solution for 7' into (128) one 
obtains an algebraic equation for 7: 



^=i«7T^r e5(4u+ "* >+16 *i- (133) 

Using this expression in the solution for B' obtained from (130) and (131), one arrives at the 
first-order equation for B 

l = 24TW (24 *'- 9 "*' 2 + 2 °"*' r <134) 

This equation determines B up to an overall multiplicative constant. This is of course a result 
of the aforementioned scaling symmetry in t, which allows us to set B hdry = 1. 

Combining these results with the solution for <j>" obtained from (128)-(129) we can isolate 
a third-order equation for the dilaton alone: 

256c//</>" - 16c// 3 (7u(j)' + 32) rht 

— ¥4, I" 



u a 2 e^ (ucf,' + 4) + 16c// u ( 5u( t>' + 12 ) ( U( t>" + V) 
x [l3uV 4 + 8u (llu 2 (p" 2 - 60c/>" - I2u<p"') + u 2 c// 3 (l3u 2 c/>" + 96) 

+2uc// 2 (-5u 3 c//" + 53n 2 0" + 36) + c// (30u 4 c// /2 - 64u 3 c//" - 288 + 32u 2 <P") 

'(135) 

At this point, we observe that the anisotropy parameter a and the undifferentiated dilaton 
only enter in eqn. (135) through the combination a 2 e2^, so that we can shift the dilaton 

0-^ = 0+^ logo, (136) 

and eliminate a from (135) altogether. The only place where a will appear explicitly is as an 
overall factor of a 2 / 7 in the eqn. (133) for 7. We also note that, since we fixed c/> bdry = 0, we 
will have 

a = el^y . (137) 
Solutions with different c/> bdry correspond then to systems with different anisotropies. 
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We can solve the equations of motion above analytically in the limiting cases of high and 
low temperature (see Appendices D and E), but for intermediate regimes we have to resort 
to numerical integration. In the next section we give some details about the procedure we 
have followed. 

A. 2 Numerics 

We have solved eqn. (135) numerically (after having shifted the dilaton as in (136)), integra- 
ting from the horizon u = u H to the boundary u = 0. 15 We expand the dilaton as 



l K + Y,^n{u-u H ) n , (138) 



ri>\ 

and similarly for the other fields, with the only restriction that J-" H = 0. Inserting the 
expansion (138) in (135) determines all the dilaton coefficients in terms of </> H and uh • These 
in turn determine the coefficients of J- and B through eqns. (133)-(134). Alternatively, one 
can substitute all the expansions into eqs. (124)-(129), and then solve at the same time for 
all the coefficients order by order. Of course, the two procedures yield the same result. Here 
we just report the expressions for a few of the coefficients: 



4 e |0H (_98304 + 34816e^ H u| + 52864e 7 ^ H u| - 336e^ H i4 + 3e 14< ^%| 



(139) 

In the following we will also need the explicit expression for J^i: 

2 1 2 



?i = — 16 + <e5^ . (140) 

4n H V / 

After integrating (135) numerically we use <f> to obtain T and B through (133)-(134). 



A. 3 Parameters 

We have seen that the solution is determined by two parameters, 4> H and n H . This was to be 
expected, since these determine the two physical parameters that the solution must depend 
on: the temperature and the anisotropy. As far as we have been able to verify with our 
numerical results, the map 

T = T(4>n, Un) , a = a(4>n, %) (141) 

is one-to-one. 

1 In practice, one actually integrates from «h — e to e, where e< 1. 
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(a) (b) 

Figure 13: (a) The metric functions for a/T ~ 4.4, with u H = 1 and H = 0. (b) The metric functions 
for a/T ~ 86, with u H = 1 and H = 0.275. 



The parameter a can be easily extracted from the boundary value of the rescaled dilaton 
bdry through (137). To compute the temperature, we have already noted that the Euclidean 
continuation of the metric (8) in the (t E , ti)-plane near u H takes the form (14). In terms of 
new coordinates 



Tx\[B^ ± 2e" 



■&= \ t E , p = y=^T^, (142) 

this piece of the metric becomes just ds E « dp 2 + p 2 d$ 2 . Regularity at u = u H then requires 
that ft have period 5$ = 2ir, and therefore that St E = Air / J 7 i\/~B^. We identify this with the 
inverse temperature: 

— tt, (l6 + 4e3M, (143) 



5t E 16iru H 

where in the last step we used (137) and (140). 

By repeating the integration procedure for different initial values of 4> H and u H , we can scan 
the (T, a) plane and construct the map (141). Roughly, we find that a depends strongly on 
cj>n, whereas T depends strongly on u H . To illustrate what the solutions look like, we present 
in Fig. 13 the plots corresponding to two different initial conditions: (cf> H = 0, n H = 1), which 
gives a moderate anisotropy-to-temperature ratio a/T ~ 4.4, and (c/> H = 0.275, n H = 1), 
which gives instead a/T ~ 86. In both cases we have checked that the metric component gtt 
in the Einstein frame is monotonically increasing towards the boundary. 



B The IR Lifshitz solution 

In this paper we have studied the finite-temperature generalization [9] of the RG flow found 
in [73] between an AdSs geometry in the UV and a Lifshitz-like anisotropic geometry in 
the IR. Neither the authors of Ref. [73] nor we have been able to write down the solution 
describing the entire RG flow analytically. However, the geometry at the IR endpoint of the 
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flow, including its finite-temperature version, was found analytically in [73]. In this appendix 
we review this solution and compare the parametrization of [73] with our Ansatz (119). 

From Sections 2.3 and 2.5 of [73] we see that the metric in string frame for the finite- 
temperature Lifshitz-like fixed point reads (with some obvious relabelings) 



ds 2 = R 2 



7 5 dv 

ra (-f(r)dt 2 + dx 2 + dy 2 ) + r^dz 2 + — 



rs/(r) 



+ R 2 rsdn 2 s5 , (144) 



where the subscript 's' stands for string frame. The radii, the blackening factor and the 
dilaton are given by 

R 2 S = ^R 2 S , f(r) = 1 - ( 11/3 , e*« = r 2 /V° , (145) 



12 s ' ' v ' \r 

where r H and <po are arbitrary constants. The forms are as in eqs. (120)-(121). As anticipated 
above, if r H = then the metric (144) in the Einstein frame is invariant under the anisotropic 
scaling 

{t,x,y} -> k{t,x,y} , z -)■ k 2/3 z , r^k~ x r. (146) 

Note however that the dilaton and the forms are not invariant. 

The coordinates inside the square brackets in the metric (144) are dimensionless. Com- 
parison with our parameterization is facilitated by making them dimensionful through the 
rescalings 

{t,x,y\^R- 1 R- 1 l % {t,x,y\, z^R^RJ^z, r^Rj 6 R 5 s r. (147) 
These bring (144) to the form 

ds 2 = (^ 7l \-f{r)dt 2 + dx 2 + dy 2 ) + ^ (148) 

where L s = (12/11) 3/5 R s and we have also rescaled r H — > (11/12) 3 r H R s , so that the black- 
ening factor still takes the form given in (145). 

The relation between r and the radial coordinate u that we have used throughout this 
paper is given by 

It is immediate to see that in terms of u the metric (148) takes precisely the form (119): 16 

L 2 / du 2 \ 
ds 2 s = -± ^ - TBdt 2 + dx 2 + dy 2 +U dz 2 + —J + L 2 s ZdVt% , (150) 

with 



3 In (119) we have set L s — 1. 
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L,\ w ' 7 /UV' 5 /L^ 2 ' 7 



and a dilaton given by 



e*(") 2/ 7^V /7 - (152) 



12/ V it 



By computing the area of the horizon, it was shown in [73] that the entropy density 

associated with the metric (148) is proportional to jV c 2 T 8//3 . By dimensional analysis, this 
implies that 

s = c ent 7V c 2 a 1 /3 T 8/3 ; ( 153 ) 

where c cnt is a constant. Note however that the constant c ont appearing above and the constant 
Cent appearing in (19) need not be the same, as emphasized in [195]. The reason is that this 
constants can be rescaled by a rescaling of the x, y, z coordinates. In the case that the solution 
approaches AdSs asymptotically, the normalization of these coordinates is fixed by requiring 
that the gauge theory metric at the boundary be canonically normalized. This normalization 
need not coincide with the normalization chosen for these coordinates in the case in which 
the solution is asymptotically Lifshitz, as is the IR geometry described by the metric (144). 

If r H = then the scalar curvature of the string- frame metric (144) diverges as r -1 / 3 
as r — > (and the same is true for the five-dimensional metric inside the square brackets). 
In contrast, all curvature invariants of the corresponding Einstein- frame metric are finite 
everywhere. Nevertheless, as pointed out in [121] and elaborated upon in [122], the hyper- 
surface r = is still problematic because tidal forces become infinite as this hypersurface is 
approached. Fortunately, if r H > then this hypersurface is hidden behind the horizon and 
the geometry is regular everywhere on and outside the horizon. 



C Relation between free energies and pressures 

In this section it will be useful to work with total energies and charges, as opposed to the 
densities used in the main text. For this purpose we consider a box of sides L x ,L y and L z 
and denote the total quantities with a bar. Setting V = L x L y L z we have: 

-v- °- G v -I <-) 

Note that, as explained in Section 4, a has dimensions of 1/length because it measures 
the number of D7-branes per unit length in the z-direction, whereas <3? has dimensions of 
energy /length 2 because it measures the energy cost per unit xy-area of introducing an addi- 
tional D7-brane extending in the xy-directions. 

The first law in terms of the total charges and allowing for independent changes in the 
lengths of the box takes the form 

dE = Tds — L x L y P z dL z — L y L z P x dL x — L z L x P y dL y + $ dd , (155) 
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from which one reads the relations 



nn (9E\ r r n ( dE\ - (BE 

V ° b J L x ,L y ,L Z ,a \ U ^z/ s,L x ,L y ,a \ . 

(156) 

Now we apply the usual extensivity argument to E = E(s, L x , L y , L z , a). For this we first 
recall Euler's theorem. This states that if a function of multiple variables / = f(cu\, ■ ■ ■ ,w n ) 
obeys the scaling law f(k ai uj\, . . . , k an u n ) = k p f(uji, . . . , w n ), then the function and its partial 
derivatives are related through 

pf(ux, ...,oj n ) = a 1 (^-) coi H h oi n ( -p-] w n • (157) 



In order to apply this to E, we note that this obeys the three scaling laws 

E(ks,kL x , Ly, L z ,a) = kE(s, L x , L y , L z ,a) , 
E(ks, L X: kLy, L z ,a) = kE(s, L x , L y , L z ,a) , 
E(ks,L x ,L y ,kL z ,ka) = kE(s, L x , L y , L z ,a) . (158) 

The first relation can be understood by imagining taking k copies of the original box and 
'piling them up' along the x-direction to make a bigger box. This scales L x — > kL x but leaves 
L y ,L z invariant. The total energy and the total entropy of course also get multiplied by k. 
However, a remains invariant because piling up the smaller boxes to make the bigger one 
along the x-direction makes the D7-branes longer, but does not change their number. An 
analogous argument with x replaced by y leads to the second relation. The third relation 
can be understood by imagining taking k copies of the original box and 'piling them up' now 
along the z-direction. The crucial difference is that in this case the number of D7-branes 
also gets rescaled as a — > ka. Any other scaling relation obeyed by E can be obtained by 
composing the three scalings above. 

Applying Euler's theorem, the relations (158) imply 



E — Ts — L x LyL z P x , 

E = TS — L X LyL Z Py , 

E = Ts - L x L y L z P z + $a. (159) 



It follows that 

P x = P y = P xy , (P z - P xy )V = §a . (160) 

Using (154) to translate to local densities and the definitions (62) and (67) finally leads to 
the relations (69)- (70). 



D High- temperature analysis 

For large values of the temperature, T 3> a, fj,, it is possible to find analytic expressions for 
the metric and the dilaton. To do so, we expand the fields (up to fourth order in a) around 
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the black D3-brane solution (17): 

u 4 

T{u) = 1- -r + a 2 T2{u) + a 4 T4{u) + 0{a 6 ), (161) 

B(u) = l + a 2 S 2 (u) + a 4 S 4 (n) + 0(a 6 ), (162) 

<t>{u) = a 2 4> 2 (u) + a A fa(u) + 0(a 6 ), (163) 

Note that only even powers can appear because of the symmetry z — > — z. We then substitute 
these expansions into Einstein's equations and solve them order by order in a, with the 
boundary conditions that all the undetermined functions vanish at the boundary u = (so 
that we preserve the AdS asymptotics and the condition 4> hdiy = 0), and that, in addition, all 
the J- n vanish at the horizon u = u H . Alternatively, one could simply expand the dilaton as 
above, solve the third-order dilaton equation (13), and then use eqs. (11) and (12) to obtain 
J- and B. However, the former method allows for an easier implementation of the boundary 
conditions. 



D.l Leading order 



We start by solving for the functions at second order, which turn out to be very simple: 



P2{u) 

B 2 {u) 
02 M 



1 



24*4 



8u 2 (ul -u 2 )- lOn 4 log 2 + (3t4 + 7u 4 ) log ( 1 + 



u 



117 



24 



10n 2 



?log 1 + 



2 ,o M 1 + ^ 



U 



(164) 



Evaluating these expressions at the horizon, it is immediate to obtain the temperature of the 
system to order a 2 : 



WBh_ 1 + « H (51og2-2) o2 + Q(fl4) 



(165) 



4-7T 7TU H 48tT 

Inverting this relation we can express the position of the horizon as a function of the tem- 
perature: 



1 51og2-2 2 4^ 

^f + ^rr a +0{a) 



(166) 



This in turn allows us to express the thermodynamical quantities in terms of the physical 
parameters a and T, rather than in terms of a and w H . Using (166) in (16) we see that the 
entropy density is given by 



JV 2 e~t*H tt 2 N 2 T 3 N 2 T 



2iru-k 



+ 



16 



0{a 4 



(167) 



The energy and pressures can also be straightforwardly extracted from the asymptotic ex- 
pansions of the fields. Using the results (35) one obtains 

A^ 2 / 3 _ 23 \ „, 4n 3vr 2 A 2 T 4 N 2 T 2 



+ 



32 



-a 2 + (9(a 4 ) , 
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i xy 
P- 



n 2 N 2 T 4 N 2 T 2 
£ + 



32 

tt 2 N 2 T 4 N 2 T 2 



32 



a 2 + 0{a 4 
a 2 + 0(a 4 



(168) 



D.2 Next-to-leading order 



We now extend the previous analysis to the 0(a) terms in (163). This is important because 
this is precisely the order at which the conformal anomaly (38) enters. 

The full expressions for JF±(u), B^u), and 04 (u) are quite cumbersome and not particu- 
larly illuminating, so we do not report them here. We limit ourselves to writing down the 
expressions for the derivatives of 04 (it) and 4>i(u): 



B' 4 (u) 



UU-r; 



1152 (u 2 - u 2 ) 2 (u 2 + ul 



2 V* 



4w 8 (371og 2 - 19) + 4it 6 it|(179 + 64 log 2) 

-4u 4 n|(301 + 807 log 2) + 3u 2 ul(188 - 952 log 2) + 64u 2 (u 2 + u|) (u 4 - 21u|) logu 2 
-4((u 2 - ul) 2 (u 2 + ul) (u 2 + llu|) log«| - 7u 2 u| (9u 4 - 10u 2 u| + 21m!) log (2« 2 ) 
-3 (54u 6 u| + 73n 4 n| + 56u 2 ul + 77«|) log (l + u 2 /u 2 ) 



+ (I5u 8 + 69n 6 n| - 394'U 4 u 4 I - 179u 2 u| - 11m|) log (u 2 + ul) 

4 

— -f- 16u 2 (u 2 + u|) log u 2 - (l6u 4 + 9u 2 ul) log (u 2 + u|) 

48(u" : — ug) (it 2 + uj) l 

+ (u 4 (3 + 37 log 2) - 7u 2 ul log (2«|) + u 2 u|(34 log 2 - 3)) 
- (I8u 2 n| + 21-u 4 ) log (1 + u 2 /ul) 



(169) 



where some constants of integrations have been fixed by requiring regularity at the horizon, 
and to writing down the equation for j-4, which can be recast in the form 



d 



J r i(x)/x 4 



dx 



288x 5 (x 2 - 1) (x 2 + if 



x 10 (19 - 189 log 2) - x 8 (29 + 189 log 2) 



+x b (13 - 211 log 2 - 48 log 2 2) + x 4 (9 - 81 log 2-48 log 2 2) 

-x 2 (12 - 30 log 2-48 log 2 2) + 48 log 2 2 - 32 (x 2 + l) (7x 4 + 3) x 4 log : 



+ (x 2 - 1) (x 2 + 1) (24Li 2 x 4 - 96Li 2 



x 2 + 1 



+ 87T 2 + 961og(2x) log (l 



-961og(l - x 2 ) log(l + x 2 ) - 9 log 2 (l + x 2 ) 
+ (63x 10 + 333x 8 + (415 + 126 log 2)x 6 - (75 - 126 log 2)x 4 

-(108 + 1261og2)x 2 + (12 - 126 log 2)) log (x 2 + l) ,(170) 

with x = u/un- 
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We also write down the values taken by all these functions at the horizon 
^4(m h ) = 0, 



4>a («h) 



369 - 8vr 2 + 912 log 2 + 354(log 2) 
6912 

3-2vr 2 -61og2 + 87(log2) 
288 



(171) 



and their asymptotic expansion near the boundary 
1 



S 4 («) 



3456 
1 



u 



915 - 40vr^ + 1611 log 2 + 1440(log 2Y + 2016 log — u 4 + 0(u b ) , 



u 



. 551 - 567 log 2 - 672 log — u 4 + 0(u b 
1152 \ u H / 

-!- (32 - 27 log 2 - 32 log —) u 4 + 0(u 6 ) . 
192 V u H J 



(172) 



Following the same logic as in the previous section, we can compute the temperature at 0(a 4 ) 
and solve for u H obtaining 



1 51og2 - 2 2 180 + 40vr 2 -121og2 - 273(log2) 



+ 



-<r + 



H ttT ' 48vr 3 T 3 " ' 13824vr 5 T 5 
The entropy (167) receives now an extra contribution 



-o 4 + 0(a 6 ) . (173) 



tt 2 N 2 T 3 N 2 T , 

+ — ~ — CL 



N 2 



,a 4 + 0(a 6 ), 



2 1 16 ~ 48tt 2 T 
while the asymptotic coefficients J-4 and B4 defined right below eqn. (31) become 
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(175) 



where we have introduced the reference scale fi required by dimensional analysis and the 
presence of the anomaly — see the discussion in Section 3. 

Substituting these expressions in (35) we find that the energy and pressures are 



37r 2 iV 2 T 4 N 2 T 2 a 2 N 2 a 4 
E = ^ + — + 
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32 ' 1536vr 2 
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8 32 ' 512vr 2 [ ~ \ a 

Eqs. (174) and (176) are precisely the results quoted at the beginning of Section 5. 



+ 0(a b ), 
+ 0(a 6 ), 
+ 0(a 6 ) .(176) 
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D.3 Check of thermodynamical relations 

First of all, we see that (176) reproduces the conformal anomaly 

iV 2 * 



-E + 2P xy + P t 



48vr 2 



(177) 



regardless of the choice of scheme, i.e. irrespectively of the values of /i and c s , 
The free energy F = E — Ts takes the form 



tt 2 N 2 T 4 N 2 T 2 a 2 



2„4 



32 



+ 



N 2 a 
1536vr 2 



8c ach -9 + 32 log 



2vrT 



+ 32 log 



+ 0(a 6 ) 



(178) 

which clearly satisfies (65) and (69). From the free energy we compute the chemical potential 
through (66): 



N 2 T 2 a + N 2 a 3 



16 



384vr 2 



8c sch -41 + 321og 



2vrT 



+ 32 log 



+ <9(a 5 ) . (179) 



and with this it is straightforward to verify that (71) holds. 

In summary, we see that all the thermodynamic relations are satisfied regardless of the 
choice of scheme despite the fact that the free energy, the pressures and the chemical potential 
are scheme-dependent. 



D.4 Scheme choice 

We conclude this Appendix by showing how the scheme choice adopted in the main text, 
c SC h = — 1) follows from the requirement that the definitions (60) and (66) agree with one 
another. 

We start from the field strength dC{u) (cf. eqn. (57)) 

d u C t xy{u) = a — -3 — . (180) 

Using the results above this equation can be integrated to next-to-leading order in a. We fix 
the integration constant so that Ct X y(uu) = and use (173) to express the result in terms of 
T and a. Expanding the result near the boundary and transforming to the FG coordinate v 
through (31), we obtain the asymptotic fall-off (58) with 

C = 1 (24avr 2 r 2 + 17a 3 ) - \c? log - \a* log Q + 0(a 5 ) . (181) 

Requiring that substitution of this expression in eqn. (60) result exactly in (179) fixes both 
the normalization chosen in (60) and c sch = — 1. 



E Low-temperature analysis 

In this appendix we consider the limit T<Co, which is also amenable to analytical treatment. 
The key observation is that at T = we must have —gtt = Qxx = 9yy and therefore FB = 1, 
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since in this limit the Lorentz symmetry in the t, x, y directions must be recovered. Either 
from the expansion of the metric coefficients (29) or from the functions themselves (28), we 
see that this implies a relation between the undetermined coefficients: 

121 

B i (a,0)+F i (a,0) = — a 4 . (182) 
57b 

Using this to eliminate B^a, 0), eqs. (35) reduce in this limit to 

P,M) = M ^r *V + |^4(..,0), (183) 
from which it follows immediately that 

E(a,0) = -P xy {a,0), P z {a, 0) = 3E(a, 0) + . (184) 

Looking at the region near the origin in Fig. 4 one can verify that these relations are satisfied 
by our numerical solution. Note also that the first relation is consistent with the fact that 
s = in this limit, since it implies that F = E. Thus we can compute the chemical potential 
as usual by differentiating with respect to a either of these quantities: 

dF(a,0)\ _ 28c sch - 1 Ar2 ^ 3 , N? dT 4 (a,0) 

T 

Requiring that this $(a, 0) also satisfies the relation P z (a,0) — P xy (a,0) = $(a, 0)a implies 
the differential equation 

a a^o)_ 7 4 = 0] 

oa 12 
which can be easily solved up to an integration constant 

^4(a, 0) = ^a 4 log Q + c int a 4 , (187) 

where as usual we have introduced the necessary scale fx. Substituting this into (183) one 
recovers the T-independent parts of (89)-(91). 

So far we have strictly set T = 0, but we can also compute analytically the first correction 
in T/a. We know in fact that when a > T our solution reproduces the Lifshitz geometry 
found in [73], and therefore that the entropy density must scale as (see the discussion below 
(153)) 

s(a, TCa) = c cnt iV c V /3 T 8 / 3 + • • • , (188) 

where c cnt ~ 3.2 is a constant that can be read off from the plot of the entropy in Fig. 2 and 
the ellipsis refers to terms of higher order in T/a. Integrating this expression with respect to 
T yields the free energy 

F(o,T«o) = -^o 1 ' 3 r 8 ' 3 + --.. (189) 
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From this and (188) one can easily obtain all other thermodynamical quantities at leading 
order in T/a. This completes the derivation of (89)-(91). 
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